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ABSTRACT 



The goal of this investigation is to design minimuin 
weight arch structures which span the distance between two 
points in two-dimensional space. An arch of unknown shape and 
variable cross-sectional width is modeled as a series of 
straight bar-beam elements. Finite Element Methods are used 
to compute the stresses in each element. Automated Design 
Synthesis (ADS) software is then used to vary the slope of 
each element and the cross-sectional width to prevent the 
yield stress of the material from being exceeded as ADS 
minimizes the arch volume to arrive at the minimum weight 
structure. Results are presented for a number of different 
loadings and boundary conditions. 
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I. 



INTRODUCTION 



A. MOTIVATION 

Since the beginning of recorded history, the arch has been 
used in the design of a wide range of structures. The arch 
was initially developed as a more reliable and efficient way 
of building a portal in a wall comprised of bricks or stone 
blocks. Over the years, it has been applied to an ever 
increasing number of load bearing structures. A significant 
reason for its popularity with architects and engineers is the 
aesthetically pleasing form of the arch. Many designers still 
enjoy the idea of an ordinary person being moved to wonder how 
a structure manages to stand with so little visible support. 

The development of the digital computer in the latter half 
of this century placed an incredibly powerful tool in the 
hands of designers who previously had to rely mainly on 
experience to select the best arch design for a given 
application. Using the high speed calculation capability of 
the computer, a designer could now use the techniques of 
numerical analysis to approximate the solution of a complex 
differential equation with a high degree of accuracy. These 
techniques could be used to model a structure under loading 
and obtain the optimum design. The arch has been the subject 
of several such studies. 
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In 1976, Farshad [ref. 1] derived two governing nonlinear 
partial differential equations for statically determinate, 
hinged-hinged arches. The system's total potential energy, 
modified by several objective functions through the use of 
Lagrange multipliers, was minimized in regards to design and 
state variables to achieve optimality and equilibrium. For an 
arch with specified span and loading, the nonlinear systems of 
equations for optimal thrust, minimum length, and minimum 
volume were presented, but not solved. 

In 1980, Rozvany et al. [ref. 2] used the Prager-Shield 
criteria to determine the optimal shape of hinged-hinge 
frames. His frames consisted of two rigidly connected 
inclined beams with a point load applied at mid-span. He 
concluded that, for a single load condition, the optimal 
structure developed either bending only or axial forces only 
in the entire arch. When the ratio 4L/D, where L is the 
length of the span and D is the constant depth of the cross- 
section, is greater than eight, only axial forces develop. In 
this case, the optimal height of the arch is L/2. When 4L/D 
is less than eight, only bending develops and the optimal 
structure is a straight beeim. The width of the cross-sections 
varied linearly from support to axis of symmetry. The author 
also found that although for a single system of point loads, 
the optimal arch consisted of straight segments, curved 
segments occurred if several alternative systems of point 
loads were considered. 
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Also in 1980, Lipson et al. [ref. 3] conducted a numerical 
study using an automated design routine to determine an 
optimal arch design with the arch shape and cross-sectional 
dimensions allowed to vary. He modelled a uniformly loaded 
parabolic arch as a system of straight segments with thin- 
walled, rectangular cross-sections. When he studied arches of 
constant depth and width, the wall thickness determined the 
optimal shape for minimum weight. The resulting shape was 
found to be a parabolic curve with a height of 0.342 times the 
span length. 

In 1988, Ang et al. [ref. 4] investigated the optimal shape 
of plastically designed non-funicular arches under a uniformly 
distributed load. The problem was solved by parameterizing 
the unspecified arch axis using spline functions and employing 
a smoothing function to approximate the non-smooth objective 
function (arch weight) . The arch considered had a rectangular 
cross-section of variable width and constant depth. It 
carried a uniform load and its supports were either hinged- 
hinged, hinged-clamped, or clamped-c lamped. The optimum shape 
of the arch was found to be a parabola with a height of 0.433 
times the span length, which appears to disagree with Lipson 's 
results . 

In 1990, Charles Scott McDavid of the Naval Postgraduate 
School [ref. 5] investigated the optimization of circular 
arches subjected to various loading and boundary conditions. 
He modelled the arches as systems of straight segments with 
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constant depth and variable width. He concluded that the 
bar/becon model was a viable technique in the approximation of 
arch structures and that the more statically indeterminate 
an arch structure was, the more efficient it was under 
identical loading. Lieutenant McDavid also proposed several 
further topics of research in this areas. One of the 
recommended topics was taken up Margaret Anne Menzies in 1991 
[ref. 6]. She investigated circular arches with varying depth 
and width. She also validated the bar /beam model as an 
approximation of an arch, but was limited by computer 
restrictions to a relatively small number of straight 
segments. This investigation looks into another area 
suggested by Lieutenant McDavid by allowing the radius of 
curvature (shape of centroidal axis) to be one of the design 
variables . 

B . PROBLEM STATEMENT 

The arch considered by this investigation has a rectangular 
cross-section with constant depth and varying thickness. The 
radius of curvature is allowed to vary to obtain the optimum 
shape of the arch. The distance spanned in the horizontal and 
vertical directions will be specified. The cross-sectional 
dimensions are assumed to be small with respect to the radius 
of curvature, which implies that the centroidal axis and the 
neutral axis coincide. 
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The governing equations for the behavior of the arch are 



the becun equilibrium equation: 

(Elv")" ^ Py(s) (1) 

and the bar equilibrium equation: 

{AEu')" = -P^{3) (2) 

The prime superscript indicates a derivative with respect to 
the independent variable's. The variables in the equations 
are defined as: 

E = Young's Modulus 

I = Cross-sectional Moment of Inertia 
A = Cross-sectional Area 
V = Lateral displacement 
u = Axial displacement 
Py = Lateral loading 
= Axial loading 
s = independent variable 

The arch is approximated by a system of straight elements. 
The local displacements are used to generate the internal 
pseudo stresses by applying the virtual load techniques 
described by Ding and Esping [ref. 7]. When the stress 
distribution is known, the volume of the arch is minimized to 
obtain a shape that keeps the developed stresses below the 
maximum allowable stress. 
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The purpose of this study is to minimize the total weight 
of an arch spanning a specified horizontal and vertical 
distance under a variety of loading conditions. The material 
of the arch is assumed to be linearly elastic, homogeneous, 
and isotropic. The slope of the arch and the cross-sectional 
width at each node are allowed to vary to obtain the optimum 
(minimum weight) structure. Automated Design Synthesis 
software [ref. 8] is used to carry out the optimization. 
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II. PROBLEM FORMULATION 



A. PROBLEM STATEMENT AND ASSUMPTIONS 

As stated previously, the goal of this investigation is to 
design a minimum weight arch which spans a pre-defined 
distance in two-dimensional space. The radius of curvature 
and cross-sectional width of the arch are allowed to vary to 
obtain the optimum shape. The following assumptions are made 
to simplify the problem: 

• The arch material is homogeneous, isotropic, and linearly 
elastic. This implies that the minimum volume arch is 
also the minimum weight arch. 

• The arch is modelled as a series of straight bar-beam 
elements whose behavior is governed by equations 1 and 2 . 

• Failure is assumed to occur when the yield strength of the 
material is exceeded. 

B . MATHEMATICAL MODEL 

The arch spans a distance in two-dimensional space from 
point A to point B (see fig. 2.1). The number of elements, 
NEL, the cross-sectional height, h, and the horizontal 
distance spanned by each element, DX, are specified. DX is 
constant and is defined as: 

DX ^ L/NEL (3) 

where L is the distance between A and B in the horizontal 
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direction. 




The design variables in this case are the cross-sectional 
base, bi, at each node and the slopes, SLPi, of the first NEL 
- 1 elements. SLPi is defined as: 

SLPi = = (Kj,, - Yi)/DX 



where Yi+j and Yi are the vertical coordinates of the nodal 
points at the ends of the i'^*' element. Since the structure 
must rise a distance H, the slope of the last element is not 
a design variable, but is fixed by whatever values the slopes 
of the previous NEL-1 elements take. 

The base dimension from one element to an adjacent element 
maintains smooth piece-wise continuity, as does a plot of the 
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spatial coordinates of the nodal points. The resulting 
elemental shape is a three-dimensional trapezoid (fig. 2.2) 
whose volume is determined by multiplying the average base and 
the height with the length of the element. The elemental 
length, 1^, is defined as: 

( 5 ) 

Therefore, the volume of the i^** element is defined as: 

Volume i ^ * h * 

where: 

+ bi)/2 

bi = nodal base dimension (6) 

h = nodal height dimension 
ij = elemental length 
i = i^ element 
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The volume of the entire arch is obtained by summing the 
individual elemental volumes. 

C. OPTIMIZATION PROBLEM 

The goal of this investigation is to obtain a minimum 
weight (volume) arch which maintains a stress distribution 
which does not exceed the yield strength of the material. The 
volume of the arch is defined as: 

im. KBL 

Volume = V Volume^ = * b * (7) 

This equation can be rewritten as: 

UEL 

Volume = h * DX * * \jl + SLPj^) (8) 

J-1 

The objective function for the optimization problem is defined 
as : 

HEL 

OBJ = ^ X * ( SLP^) ^ (9) 

Since h and DX are constant, minimizing the objective function 
will minimize the arch volume. 

Constraints are imposed upon the objective function and 
design variables to ensure the applicability of the governing 
equations and to prevent failure by yielding. The strength 
constraint can be simply stated as follows: 
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( 10 ) 



or in normalized form-. 

- 1.0 s 0.0 

vfheiex 

aj = maximum stress at i^ node 
Sy = material yield strength 

The stress distribution in the arch is obtained by Finite 
Element Methods. This process is described in Chapter IV. 

The constraints placed on the base dimensions of the arch 
are used to ensure the validity of the bar and beam 
equilibrium equations (equations (1) and (2)). Limiting the 
base relative to the constant height prevents the structure 
from behaving like a shell ar a deep curved beam. The 

following limits are imposed through the use of the ADS side 
constraints on b^ ; 

Q.l*h ^ ^ 3*h (11) 

It is also important to choose a cross- sectional height 
dimension that is small relative to the length of the entire 
structure. 

One limit on the cross-sectional dimensions that was not 
considered by this investigation is that imposed by the 
requirement for elastic stability. The obvious importance of 
avoiding failure by buckling is an interesting area for future 
research. 
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III. OPTIMIZATION ANALYSIS 



The computer optimization is performed by ADS, a general 
purpose numerical optimization program containing a wide 
variety of algorithms. The program minimizes one function, 
the objective, subject to bounds imposed on the design 
variables by constraint functions. The solution of the 
general problem is separated into three basic levels which are 
strategy, optimizer, and one-dimensional search. For this 
constrained minimization problem, the ADS user's guide [ref. 
8] recommends a Sequential Linear Programming strategy, a 
Modified Method of Feasible Directions optimizer, and a Golden 
Section Method one-dimensional search. 

A. STRATEGY 

Sequential Linear Programming (SLP), linearizes nonlinear 
objective and constraint functions, then obtains a solution to 
the approximation of the problem by linear progrcunming methods 
(see fig. 3.1). The problem is linearized again about the 
design point yielded by the solution to the previous 
approximation and resolved. This process is repeated until a 
precise solution is obtained [ref. 9]. 

The nonlinear functions are linearized through the use of 
a first-order Taylor Series expansion in the following manner; 
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Minimize: F{X) - F{X^) + VF(-Y^) ♦ hX 
Subject to: gjiXf,^ + Vgj{x^) * bx ^ o 
where: bx = X - X^^ 

j = constraint 

The zero subscript identifies the point about which the Taylor 
series expansion is performed. At the initial design, the 
objective and constraints are linearized to give straight line 
representations of the functions. 

In an under-constrained problem like this one, where there 
are fewer active constraints than design variables, this 
method sometimes performs poorly. This occurs because the 
linear approximation may be unbounded. This problem is dealt 
with by imposing move limits on the linear approximation as 
shown in figure 3.1. This ensures that the optimum is 
eventually reached within a tolerance of the move limits. In 
practice, these move limits are reduced during the design 
process so that a solution is found with the desired accuracy. 
SLP tends to converge rapidly to a solution, but while the 
solution of the linear problem is near the nonlinear optimum, 
it is in the infeasible region. 



B. OPTIMIZER 

The Modified Method of Feasible Directions for constrained 
minimization is used to solve the approximate optimization 
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rZ^^\^LE, REGION 




OBJECTIVE ryNcriPN contours fp) 

Trmmn CONSTRAINT Pl/NCTION 

UH6AR. approximation 

Figure 3.1: Sequential Linear Progrcunming 



sub-problem (eg. 12). In this method, a search direction 
vector is first found. The vector containing the design 
variables is then updated by moving in the direction of the 
search vector in the following manner (refer to fig. 3.2): 
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( 13 ) 



X = X . + a*S 

"^q-l 

The scalar a* defines the distance moved in the direction of 
the search vector in design variable space and q represents 
the iteration number. From the initial design, the objective 
function is moved in the direction of steepest descent in 
design variable space. This will continue until a constraint 
is encountered. Once this happens, a new search direction is 
obtained by solving the sub-problem: 



Maximize: * S 

Subject to: ^g^ix) *5^0 j ^ J 
S * S ^ 1 



where: 

^ s vector whose elements are the design variables 
F - the objective function 
S = the search vector 
J = number of active constraints 



(14) 



The search direction will follow the constraint, but will 



allow the design to leave the constraint boundary if the 
objective function can be further reduced. If the scalar 
product of the gradient of each critical constraint with the 
search vector is less than zero, the search vector is moving 
away from the currently active constraint. This constraint is 
then dropped from the set of active constraints. If the 
search vector is the null vector, or numerically small, the 
optimization is terminated because this indicates that the 
Kuhn-Tucker conditions for optimality have been met [ref. 9]. 
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Figure 3.2: Modified Method of Feasible Directions 

C. ONE-DIMENSIONAL SEARCH 

The scalar a* in equation (13) is found by the Golden 
Section Method. Basically, this determines how far the 
optimizer will search in the direction of the search vector. 
The one-dimensional search attempts to find the optimal value 
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of a* which will result in the minimum value of the objective 
function. It accomplishes this by progressively bracketing 
the minimum of the function by comparing the values of the 
function at smaller and smaller limits. 

D. ADS PROGRAM PARAMETERS 

Reference 8 contains a table of ADS program parameters. 
It is possible to change these parcuneters from their default 
values with the first call to ADS using an INFO value equal to 
-2. The purpose of modifying these parameters is to try to 
obtain a better optimal design by "fine tuning" the program. 
For instance, modifying the constraint tolerance, CT, can 
ensure that the active constraints will be closer to the 
desired value. However, this can also mean that it will take 
much longer for the solution to converge. The only pareuneter 
modification done in this investigation was to inhibit the 
auto-scaling function of ADS. 
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IV. STRESS ANALYSIS 



In order to accomplish the goal of designing a minimum 
weight arch, a FORTRAN computer was used to evaluate the 
objective function and constraints and to invoke the ADS 
software. The strength constraints require that the stress at 
each nodal point remain below the yield strength of the 
material. Determining the stress distribution in the arch is 
a complex problem, because of the need to be able to analyze 
a statically indeterminate structure. Finite Element Methods 
are used to calculate the displacements at each nodal point. 
The displacements are then related to the nodal forces and 
moments by the direct stiffness method. Once the forces and 
moments are known, the stresses can be evaluated with a 
knowledge of the cross-sectional dimensions. 

A. STRESS DEVELOPMENT 

For the purposes of this investigation, the arch is 
assumed to fail if the total stress at any nodal point exceeds 
the yield stress of the material. The total stress is defined 
as the sum of the normal stresses due to bending and axial 
loading. 
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wheiei 

= total normal stress (15) 

= normal stress due to bending 
Og = normal stress due to axial load 

The absolute values of the bending and axial stresses are used 
because at either the top or bottom extreme fiber the stresses 
will be additive. In other words, the yield stress is assumed 
to be the same in compression and tension. There will be 
shear stresses developed in the arch, but the geometric limits 
imposed by the side constraints ensure that these stresses 
remain negligible (see Appendix A) . 

To calculate the two normal stress components, the arch is 
modelled as a series of straight elements, which can be 
considered to behave as beams to calculate the bending 
stresses and as bars to calculate the axial stress. The 
bending stress in an element is defined as: 







where: 

M = bending mcanent 



(16) 



c = distance frcmi centroidal axis to extreme fiber 



I * Moment of Inertia of cross-section 



The moment, M, is calculated by: 
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( 17 ) 



M = Elv^^ 



where i 

E = Young ^8 Modulus 
V = lateral displacement 



By substituting equation (17) into equation (16) the following 
equation is obtained: 



In a similar manner, the axial stress is calculated. The 
stress due to axial loading is defined as: 



The axial force, F, is determined by the following equation: 

F = AEu' 

where: (20) 

u = axial displacement 

By substituting equation (20) into equation (19) f the 
following result is achieved: 



Combining equations (18) and equation (21) results in: 



= Ecv^' 



( 18 ) 



where: 

F - axial load 
A = cross-sectional area 



( 19 ) 




( 21 ) 



= E{CV^^ + U') 



( 22 ) 
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The approximate values of the lateral and axial 
displacements are obtained through the use of the Galerkin 
Finite Element Method. With the use of these values, the 
stresses at each nodal point can be calculated. 

B. FINITE ELEMENT BEAM EQUATION DEVELOPMENT 

The Galerkin Finite Element Method is used to transform 
the fourth-order differential equation governing static 
displacement of the beam into a system of linear algebraic 
equations. In order to maintain continuity of slopes and 
displacements from element to element, a family of cubic shape 
functions is introduced [ref. 5 and ref. 6]. An approximate 
solution for lateral displacement is defined as: 

V » V 

where: 

V = exact solution in continuous space (23) 

= approximate solution in discrete space 

N = vector containing shape functions 

V = vector containing lateral displacements and slopes 

The next step is to form a measure of the error of the 
approximation, or residual: 

R = [EIV^'] - Py{s) (24) 

where Py is the lateral load and s is the independent 
variable. Substituting equation (23) into equation (24) 
yields the following equation: 
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( 25 ) 



R= [EI(N^9)"]" - Pyis) 
The residual can be minimized if : 



f m{3)ds = 0 

J D 



( 26 ) 



where the quantity on the right hand side of the equation is 
the null vector. Substitution of equation (25) into equation 
(26) gives the following result: 



- jjipy(s)d3 = 0 



( 27 ) 



Integration by parts is performed twice on equation (27), 
resulting in: 



N[EI{N^X>-)'']'\b - (E)'EHN^^)"\b 
+ j "EUN^^) "ds - jNpy{s)ds = 0 



( 28 ) 



where I ^ indicates an evaluation at the boundary points of the 
structure. Since the displacement vector is a constant, 
equation (28) can be rewritten as: 



jj.Sri"Ei{N'^)"ds^ - jSpy(s)ds = 0 



( 29 ) 



From the beam equilibrium equation, the shear force, V, is 
defined as: 



V = EIv 



./// 



( 30 ) 



and the moment, M, by: 
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(31) 



M = Elv" 

Therefore, the boundary term load vectors are defined as: 

V = 

and: 

M= 

A system stiffness matrix is defined as: 

Kg • f (N) "ei{S'^) "ds 

J D 

and a system force vector by: 

Ft, = ^Spy{3)ds 

Substitution of equations (32) through (35) into equation (29) 
yields the following system of linear equations: 

V\g - M\g *Kg1^-Ft,= 0 (36) 

If a load vector of internal and external applied loads is 
defined as: 

= h* "U - (37) 

equation (36) reduces to: 

(38) 

The global bending stiffness matrix is constructed from the 
union of the elemental bending stiffness matrices and the 
global bending force vector is constructed from the union of 
the elemental bending force vectors. 



(32) 

(33) 

(34) 

(35) 
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C. FINITE ELEMENT BAR EQUATION DEVELOPMENT 

The development of the Galerkin Finite Element Method for 
the bar equilibrium equation is done in a manner similar to 
that for the beam equation. However, since the bar equation 
is a first-order differential equation, a family of linear 
shape functions can be used to maintain continuity of 
displacement from element to element [ref. 5 and ref. 6]. 
Once again, an approximate solution for the axial 
displacement, u, is formed as follows: 

u » a = N'^u 

where: 

u = exact solution in continuous space 
Q = approximate solution in discrete space 
N = vector containing shape functions 
U = vector containing axial displacements 

The residual is defined as: 

Ris) = [AEiN^U)']' * pjs) 

The residual can be minimized if: 

f N[R(s)]ds = 0 

J D 

Substituting equation (40) into equation (41) results in the 
following equation: 

^N[AE(S^Q)']'di3 * lj}pjs)ds = 0 (42) 

Unlike the beam equation development, only one integration by 
parts is performed to yield; 



(39) 



(40) 



(41) 
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A£W(W’^u)'L - f {N)'[AE{if'^a)']ds + f Np^(s)ds = 0 (43) 

J D J D 

where |b indicates an evaluation at the boundary points of the 
structure. Since the axial displacement vector is a constant. 



equation (43) can be rewritten as: 

N{AEi}^)'a\s - f (N)'[AE(N^)']dsU + f Np^ = 0 (44) 

J D J D 

From the bar equilibrium equation, the axial force, F, is 
defined as: 

F = AEu' (45) 

Therefore, the boundary term load vector vectors are defined 
as : 

p = AEN{N'^ya\g (46) 

A stiffness matrix is defined as: 

( (ih'[AE{N'^)]ds (47) 

V D 

and a system force vector by: 

F^^fNpJs)ds (48) 

V D 

Substituting equations (46) through (48) into equation (44) 
yields the following equation: 

p - Kjfi + = 0 (49) 

By defining a new load vector as: 
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( 50 ) 



equation (49) reduces to: 




( 51 ) 



The global axial stiffness matrix is constructed from the 
union of the elemental axial stiffness matrices and the global 
axial force vector is constructed from the union of the 
elemental axial force vectors. 

D. THE ELEMENTAL STIFFNESS MATRIX 

The global Galerkin FEM stiffness matrices in equations 
(38) and (51) are constructed from the union of elemental 
axial and bending stiffness matrices. A single beam element 
has four degrees of freedom, which yields a 4x4 elemental 
stiffness matrix. A single bar element has two degrees of 
freedom, which yields a 2x2 elemental stiffness matrix. These 
elements can be combined into a single bar-beam element (fig. 
4.1) which has six degrees of freedom [ref. 5 and ref. 6]. 




Figure 4.1: Bar-Beam Element - Degrees of Freedom 
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This results in a 6x6 elemental stiffness matrix. The 



elemental displacement vector can be expressed as: 

where for the i ^ element: 

= axial displacement at local node l 
^2^^ = lateral displacement at local node 1 
63-^' = beam slope at local node 1 
54-^' = axial displacement at local node 2 
= lateral displacement at local node 2 
= beam slope at local node 2 



The elemental force vector can be expressed as: 



where for the i ^ element: 

= axial force at local node l 
f2^^ = lateral force at local node 1 
f^^* = mcmient at local node 1 

= axial force at local node 2 
fg-^' = lateral force at local node 2 
= m<xnent at local node 2 



Therefore, the combination of the Galerkin bar and 
equations yields: 






( 52 ) 



( 53 ) 



beam 

( 54 ) 



where the elemental stiffness matrix takes the following form: 



AE 



0 



0 



0 



F' = 



AE 



12 



El 

El 



El 
El 

0 

- 12 -^:^ - 6 -^ 

2 / 12 



■f » 



El 

El 



, El 



0 -12 
0 - 6 - 
^ 0 



El 

El 



12 

-6 



El 

El 



0 



. El 

' V 

£:r 



( 55 ) 



It is apparent from the form of the elemental stiffness matrix 
that the bar and beam have uncoupled behavior. 



E. COORDINATE TRANSFORMATION OF THE ELEMENTAL SYSTEM OF 
EQUATIONS 

Each element in the arch will have a unique orientation 
with respect to the global x and y axes. In order to solve 
the global system of equations, it is necessary to transform 
each elemental Galerkin equation into global coordinates. The 
global reference coordinate system is defined as the 
horizontal and vertical axes of the entire arch. Each element 
makes an angle Oi with the horizontal axis and an angle J3i with 
the vertical axis. 

The local displacements and forces (indicated by prime 
superscript) are defined as follows: 
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= fii-^cos (a^) + 

^2“^^ = “6i-^cos(Pj) + 62^^06 (ttj) 
63 ^' = 63 ^ 

64-^' = 64-^cos(aj) + d5-^cos(Pj) 
^ 5 ^^ “fi4-^cos(Pj) + 65 '^cos(aj) 

and: 

= fi-^cos(aj) + f2^cos(Pi) 
±2^* = ■*■ f2^cos{a^) 

f = f i 
■^3 -^3 

- f^^cosiai) + f5-^cos(Pj) 
fg-^^ = -f4-^COS(Pj) + fg-^COS(aj) 



( 56 ) 



(57) 



Therefore a transformation matrix can be defined as: 






cos(aj) cos(Pj) 0 0 0 0 

-cos(Pj) cos(aj) 0 0 0 0 

0 0 1 0 0 0 

0 0 0 cos(aj) cos(Pj) 0 

0 0 0 -COB(Pj) COB(Ctj) 0 

0 0 0 0 0 1 



(58) 



This reduces equations (56) and (57) to: 

5i' 



(59) 
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and: 



f=F8-‘ (60) 

where : 

(f^)’-= (62) 

The transformed elemental stiffness equation becomes: 

F'r'fi' = (63) 

Multiplying both sides of equation (63) by the inverse of the 
transformation matrix, which is an orthogonal matrix, yields: 

(P)^F'(r^)5^ = (64) 

where the elemental stiffness matrix in terms of global 
coordinates is defined as: 



ki = {r^)'Tfc^'(r^) (65) 
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F. SOLUTION 



The global stiffness matrix and global force vector are 
constructed from the union of the elemental stiffness matrices 
and elemental force vectors, which yields the following 
equation: 

FA = F 

whezei 

K = global stiffness matrix ^55^ 

A = global displacement vector 
F = global force vector 

By multiplying both sides of the equation by the inverse 
of the global stiffness matrix, the global displacement vector 
can be calculated. The global displacements are then 
transformed back into the local displacements and can be used 
in conjunction with the elemental stiffness matrix to 
calculate the local forces: 

( 67 ) 

These local forces and bending moments are used to 
calculate the stress at each nodal point. The stresses of 
internal global nodal points are averaged since local nodal 
point 2 of the i^^ element is the same as local nodal point 1 
of the (i+l)^** element. 
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V. PROGRAM DESCRIPTION 



The Fortran 77 code originally used in references 5 and 6 
to analyze circular arches was modified to analyze the non- 
circular arches of interest in this study. A copy of the 
program is included in Appendix B. The program reads 
information from an input file which provides the following 
user-supplied data: 

L = horizontal distance spanned by arch 

H = vertical distance spanned by arch 

HGT = depth of cross-section 

YOUNG = Young's Modulus 

YIELD = yield strength 

NEL = number of elements 

ISTRAT = ADS parameter which designates strategy 

lOPT = ADS parameter which designates optimizer 

ZONED = ADS parameter which designates one-d search method 

IPRINT = ADS print control parameter 

IGRAD = ADS parameter which designates method of gradient 
calculation 

DV__BG = initial value of design variable 

DV_LO = lower side constraint on design variable 

DV_UP = upper side constraint on design variable 

CLAN = node at which concentrated load is applied 

FX = concentrated horizontal force (positive to the right) 
in pounds 
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FY = concentrated vertical force (positive upwards) in 
pounds 

FM = concentrated moment (positive counter-clockwise) in 
pound-inches 

FA = distributed load in pounds per inch 

BX_ = boundary condition on horizontal displacement at 1, 
first node, or 2, last node (= 1 - fixed, = 0 - free) 

BY_ = boundary condition on vertical displacement at 1, 
first node, or 2, last node (= 1 - fixed, = 0 - free) 

BM_ = boundary condition on rotation at 1, first node, or 2, 
last node (=1, - fixed, = 0 - free) 

LABEL = string character identifying particular case study 

The program then utilizes several subroutines to perform the 
FEM analysis and call the ADS program to perform the 
optimization. 

An outline of the program, ARCHOPT, and its subroutines is 
provided in fig. 5.1. The first subroutine called by the main 
program is OPTIMIZATION_TOOL, which establishes the ADS 
parameters before the first call of the ADS program. The 
first call serves to override some of the default parameters 
in order to "fine-tune" the program. After ADS is called, 
OPTIMIZATION^TOOL calls subroutine EVAL to evaluate the 
objective function and constraints, which are functions of the 
design variables. 

Subroutine EVAL calls subroutine ARCH_STRESS, which in 
turn calls subroutines FORM and FORCE_VECTOR . The latter 
subroutines are used to form the global stiffness matrix and 
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the global force vector, respectively. Subroutine BNDRY 
imposes the user-supplied boundary conditions by modifying the 
stiffness matrix and the force vector. An equation solver 
from the IMSL library, L2ARG, is used to solve the Galerkin 
FEM equation by matrix inversion. Once the global 
displacements are known, they can be used to calculate the 
nodal stresses, which is done in subroutine STRESS. Once the 
design constraints are evaluated for the initial design, the 
problem is returned to ADS, where an updated design is chosen. 
This process continues until the solution converges. 

Once the termination criteria for optimization is 
satisfied, the main program calls subroutine ARCH_OUT, which 
creates an output file, ARCHOUT, which contains the problem 
parameters, the optimized design, and the arch volume. Copies 
of the output files for all case studies are provided in 
Appendix D. 

The methodology of the program was thoroughly validated in 
references 5 and 6. Once the modifications were made, the 
program was used to analyze several statically determinate 
structures (i.e. a cantilever beam) to verify that the 
approximate values of stress and displacement came close to 
the exact values. The results confirmed that the FEM analysis 
provided an excellent approximation, with less than 0.8% error 
for a twelve element model of a cantilever beam with a 
lateral concentrated load at the free end. 



34 




Figure 5.1: Program ARCHOPT Structure 
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VI . CASE STUDIES 



The optimized designs for eleven different cases of 
loading and boundary conditions are presented. Cases (1) 
through (4) are hinged-hinged arches which span equal 
distances in the horizontal and vertical directions. Each 
case has a different loading. Cases (5) through (8) are 
arches which span equal horizontal and vertical distances, but 
which have different boundary conditions. Cases (9) and (10) 
are hinged-hinged arches which span different distances in the 
vertical direction. Plots of arch shape, base dimensions, and 
stress are provided for each case. It should be note that the 
distributed loads are applied in a direction normal to each 
element. This would closely approximate a situation where the 
arch was under hydrostatic, or pressure loading. Twelve 
elements were used to model the arch in all cases. This 
implies that there are thirteen nodal points and thirteen 
stress constraints. A Young's Modulus of 30,000,000 psi and 
a Sy of 52,000 psi are used in all cases. The output files 
from progreun ARCHOPT are presented in Appendix C. 
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A. CASE 1: HINGED-HINGED ARCH WITH DISTRIBUTED LOAD 

This design has nine active stress constraints and 
thirteen active side constraints (lower limit on base). The 
stress due to axial loading is higher than the bending stress 
at all nodes. This indicates that the minimum area cross- 
section is better able to withstand an axial load than a 
bending moment. A predominantly axial stress distribution is 
more efficient since the entire cross-section is stressed to 
the same degree. In the case of a bending stress 
distribution, the material has to be strong enough to 
withstand the extreme fiber stress, which means that the 
internal part of the cross-section will be stressed below the 
yield strength. 

• H = 32.0 inches 

• L = 32.0 inches 



• h = 1.5 inches 

• Distributed load = -100 pounds/inch 

• Volume = 10.324 cubic inches 
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Case 1: Arch Shape and Width 




0 ) 

OT 

CO 



PQ 
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Figure 6.2: Case 1 - Arch Shape and Width 
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Case 1 : Stress 




Axial Stress — ^ Bending Stress Total Stress 



Figure 6.3; Case 1 - Stress 
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B. CASE lA: HINGED-HINGED BEAM WITH DISTRIBUTED LOAD 

This case is presented for comparison to Case 1. The 
slopes of the elements are not allowed to vary, so the arch 
remains a straight beam. There are seven active stress 
constraints and no active side constraints. The bending 
stresses dominate and the axial stresses are negligible. The 
volume of this arch is more than six times greater than that 
for Case 1. 

• H = 32.0 inches 

• L = 32.0 inches 

• h = 1.5 inches 

• Distributed Load = -100 pounds/inch 

• Volume = 68.387 cubic inches 
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Case lA: Arch Shape and Width 




Y-Coord — Base 



Figure 6.5: Case lA - Arch Shape and Width 



Case 1 A: Stress 




Bending Stress 



Figure 6.6: Case lA - Stress 



42 



C. CASE 2: HINGED-HINGED ARCH WITH CONCENTRATED LOAD 



This arch has a concentrated vertical load applied 
downward at node 7. There are seven active stress constraints 
and thirteen active side constraints (lower limit on base) . 
The arch has taken the shape of two straight lines meeting at 
the center node. At eight of the thirteen nodes the axial 
stress dominates. 

• H = 32.0 inches 

• L = 32.0 inches 

• h = 1.5 inches 

• Concentrated Load = 2000 pounds (downward) 

• Volume = 10.239 cubic inches 
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Case 2: Arch Shape and Width 




Y-Coc«L ^ BtK 



Figure 6.8; Case 2 - Arch Shape and Width 



Case 2: Stress 
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Figure 6.9; Case 2 - Stress 
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D. CASE 3: HINGED-HINGED ARCH KITH CONCENTRATED LOAD 

This design has five active stress constraints and twelve 
active side constraints (lower limit on base). The arch has 
taken the shape of two straight lines meeting at node 7, where 
the concentrated load is applied. Once again, this design 
shows a preference for axial loading over bending moments. 

• H = 32.0 inches 

• L = 32.0 inches 

• h = 1.5 inches 



• Concentrated Load = 2828 pounds (applied down and to the 
right ) 

• Volume = 10.842 cubic inches 
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Case 3: Arch Shape and Width 




♦ Y-Coofd. ^ Biik 



Figure 6.11: Case 3 - Arch Shape and Width 



47 




Case 3: Stress 
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Figure 6.12: Case 3 - Stress 



48 



E. CASE 4: HINGED-HINGED ARCH KITH CONCENTRATED MOMENT 

This design has one active stress constraint and twelve 
active side constraints (lower limit on base). The arch 
assumes the shape of a straight line, which is the minimum 
shape for a feasible design. The maximum base is at node 7, 
where the concentrated moment is applied. 

• H = 32.0 inches 

• L = 32.0 inches 

• h = 1.5 inches 

• Concentrated Moment = 6000 pound-inches 

• Volume = 10.207 cubic inches 
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Case 4: Arch Shape and Width 




Y-Cooid ^ B«c 



Figure 6.14; Case 4 - Arch Shape and Width 
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Case 4: Stress 




Rcadag Strew 



Figure 6.15: Case 4 - Stress 
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F. CASE 5: ROLLER-HINGED ARCH WITH DISTRIBUTED LOAD 

This design has ten active stress constraints and twelve 
active side constraints (lower limit on base). The arch has 
assumed a nearly circular shape. In this case, the bending 
stress dominates throughout the arch. It appears that the 
arch has attained the maximum degree of axial loading that can 
be achieved with these boundary conditions. 

• H = 32.0 inches 

• L = 32.0 inches 

• h = 1.5 inches 

• Distributed Load = -100 pounds/inch 

• Volume = 11.329 cubic inches 
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Case 5: Arch Shape and Width 




Y-Coofd ^ I 



Figure 6.17: Case 5 - Arch Shape and Width 
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Base (in.) 



Case 5: Stress 
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Figure 6.18; Case 5 - Stress 
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G. CASE 6: ROLLER-FIXED ARCH WITH DISTRIBUTED LOAD 

This design has four active stress constraints and ten 
active side constraints (lower limit on base). The bending 
stresses dominated for the most part, but the axial stresses 
are significant. The maximum base is at the fixed end, as 
expected. 

• H = 32.0 inches 

• L = 32.0 inches 

• h = 1.5 inches 

• Distributed Load = -100 pounds/inch 

• Volume = 14.431 cubic inches 
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Case 6: Arch Shape and Width 




Y-Cood 



Figure 6.20: Case 6 - Arch Shape and Width 
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Case 6: Stress 
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Figure 6.21: Case 6 - Stress 
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H. CASE 7: FIXED-FREE ARCH WITH DISTRIBUTED LOAD 

Using the cross-sectional depth of the previous case 
studies resulted in a design that consistently violated the 
stress constraints at the nodes at and close to the fixed end 
of the arch. Therefore a cross-sectional depth of 3 inches 
was used. The resulting design has seven active stress 
constraints and six active side constraints (lower limit on 
base) . The arch has taken the shape of a tapered beam with a 
maximum base at the fixed end. The bending stresses dominate 
and the axial loads are negligible. Because of the lack of 
support at the free end, the arch is unable to attain axial 
loading in the elements. 

• H = 32.0 inches 

• L = 32.0 inches 

• h = 3.0 inches 



• Distributed Load = -100 pounds/inch 

• Volume = 72.614 cubic inches 
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Case 7: Arch Shape and Width 
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Figure 6.23: Case 7 - Arch Shape and Width 



0 ) 

w 

(0 

CQ 



59 



Case 7: Stress 
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Figure 6.24: Case 7 - Stress 
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I. CASE 8: FIXED-FIXED ARCH WITH DISTRIBUTED LOAD 

This design has six active stress constraints and eleven 
active side constraints. The arch curves outward slightly in 
the center. With the exception of the nodes at the ends, 
axial stresses dominate in the arch. The maximum base 
dimensions occur at the fixed ends. 

• H = 32.0 inches 

• L = 32.0 inches 

• h = 1.5 inches 

• Distributed Load = -100 pounds/inch 

• Volume = 11.308 inches 
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Case 8: Arch Shape and Width 
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Figure 6.26: Case 8 - Arch Shape and Width 



62 



Case 8: Stress 
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Figure 6.27: Case 8 - Stress 
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J. CASE 9: HINGED'HINGED ARCH WITH DISTRIBUTED LOAD 

This case differs from Case 1 only in the distance spanned 
in the vertical direction , H. In this case, H = 18.475 
inches. There are eleven active stress constraints and 
thirteen active side constraints (lower limit on base). The 
arch has assumed a curved shape similar to Case 1. Once 
again, there is a clear preference for axial loading. 

• H = 18.475 inches 

• L = 32.0 inches 

• h = 1.5 inches 

• Distributed Load = -100 pounds/inch 

• Volume = 8.391 cubic inches 
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Case 9: Arch Shape and Width 




H*- Y-Coord ^ 



Figure 6.29: Case 9 - Arch Shape and Width 
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Case 9: Stress 
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Figure 6.30; Case 9 - Stress 
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K. CASE 10: HINGED-HINGED ARCH WITH DISTRIBUTED LOAD 

This case differs from the previous case only in the 
distance spanned in the vertical direction, H. There are 
eight active stress constraints and twelve active side 
constraints. Once again, the arch curves outward and the 
axial stresses dominate. 

• H = 55.426 inches 

• L = 32.0 inches 

• h = 1.5 inches 

• Distributed Load = -100 pounds/inch 

• Volume = 15.149 cubic inches 
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Case 10: Arch Shape and Width 
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Figure 6.32: Case 10 - Arch Shape and Width 



68 



Case 10: Stress 




-•-AriilSdwi -«-Bea£ii| Sum -w- Total Smi 



Figure 6.33: Case 10 - Stress 
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VII. CONCLUSIONS 



The conclusions of this study are: 



• The objective function appears to have local optimum 
points, because the initial shape of the arch affected the 
output from ADS. This required a trial-and-error approach 
to find the true optimal design for each case. 

• In the optimum designs, the base dimensions were reduced 
to the maximum extent possible, and the arch took a 
rounded shape to decrease bending moments and increase 
axial loading. One exception to this was the case of the 
fixed-free arch, which was unable to attain axial loading 
and took the shape of a straight, tapered beam. The only 
other exception was Case (4) , the hinged-hinged arch with 
a concentrated moment. 

• The smooth, curved shapes that most of the arches assumed 
confirmed that modelling an arch as a series of straight 
elements was a valid approach as long as the number of 
elements was relatively high (greater than eight) . 

• When using the combination of strategy and optimizer 
(S.L.P. and M.M.F.D.) that this investigation employed, it 
is critical that the initial design be feasible. Since 
the initial shape in all cases was a straight beam, the 
initial base dimensions had to be large enough to allow 
the beam to be able to withstand the bending stresses that 
would dominate in such a structure. If the initial design 
was not feasible, the resulting output from ADS had many 
violated constraints. 

• The most inefficient method of spanning the distance from 
point A to point B was the cantilever beam (Case (7)). 

• Case Studies (5) through (8) confirmed the conclusions of 
references 5 and 6 that, for a given loading, the more 
statically indeterminate the structure, the greater the 
efficiency. However, Case (1) contradicted this 
conclusion since it is less statically indeterminate than 
Case (8) but has a slightly smaller volume. 
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Suggested areas for future research are; 



• Optimization of an arch with a more complex cross-section 
(i.e. thin-walled tubes, WF's, etc.) 

• Optimization of arches with geometric constraints which 
prevent global buckling and local crippling. 
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APPENDIX A: JUSTIFICATION FOR OMITTING SHEAR STRESSES 

(The following appendix is taken from ref. 5) 



The shear stress distribution through a beam of 
rectangular cross-section has a parabolic distribution along 
the depth of the member. The maximum shear stress, located at 
the neutral axis of the beam, is (from ref. 10): 




wherei 

= maximum shear stress 
V = shear force 
A = cross-sectional area 



(B.l) 



The maximum normal stress due to bending is given by: 



Me 

I 



where: 

= maximum normal stress 
M = bending moment 

c = distance frexn N.A. to extreme fiber 

hh^ 

I = -f— = mcxnent of inertia 



(B.2) 



Redefining the normal stress in terms of the cross- 
sectional dimensions yields: 
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= «(A/2)/(i)AVl2) 



(B.3) 



o n = 6M/hA 



The ratio of the maximum shear stress to the maximum 
normal stress due to bending is denoted by r and given by: 



_ *max 



(B.4) 



Substituting equations (B.l) and (B.3) into equation (B.4) 
yields; 



r = il,5V/A)/{6M/hA) 
or: 




(B.5) 



For the cases investigated in this study, the maximum value r 
can attain is when the loading is that of a uniformly 
distributed load, Py. In this case: 



V- PyL 

and: (B.6) 

M = PyL^/2 



which upon substitution into equation (B.5) yields: 



73 



T = {p^)h/i{p^^/2) 

or: 




(B.7) 



The use of the beam equation requires that length of the 
beam to be at a minimum ten times the height, or in other 
words : 

L ^ 10/3 (B.8) 

To maximize the value of r, let L equal lOh, the minimum 
allowable length. Substituting this value of L into Equation 
(B.7) yields: 

r ^ h/2 ( 10 / 3 ) 

or: (B.9) 
r ^ 1/20 

Therefore, the maximum shear stress accounts for less than 
5% of the bending stress developed in the structure. Five 
percent is high considering this analysis over-assumed the 
value of the shear stress by assigning the maximum shear 
stress to the entire cross-section of the beam. Moreover, at 
the outermost fibers where the normal stress is a maximum, the 
shear stress is zero. Therefore, under the circumstances of 
this study, the addition of shear stresses was deemed to be 
unwarranted . 
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file: archopt fortran a 



APPENDIX B 



Program ARCHOPT 



PROGRAM ARCHOPT 

Cx X 

cx ARCH OPTIMIZATION ANALYSIS CODE x 

Cx X 

Cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

c 

C ALPHA TRANSFORMATION ANGLE OF ELEMENT (ANGLE TO X-AXIS) 

C SAVE THE AVERAGE BASE DIMENSION ACROSS AN ELEMENT 

C BASE ARRAY CONTAINING THE ELEMENTAL BASE DIMENSIONS 

C BASEL ARRAY CONTAINING THE ELEMENTAL BASE DIMENSIONS LOV/ER 

C SIDE CONSTRAINT 

C BASEU. .. .ARRAY CONTAINING THE ELEMENTAL BASE DIMENSIONS UPPER 
C SIDE CONSTRAINT 

C BETA TRANSFORMATION ANGLE OF ELEMENT (ANGLE TO Y-AXIS) 

C B_1 BOUNDARY TERMS APPLIED AT END ”1" 

C B_2 BOUNDARY TERMS APPLIED AT END ”2" 

C Cl,. . ,C5. CONSTANTS RELATED TO ELEMENT STIFFNESS COEFFICIENTS 
C CLAN..... CONCENTRATED LOAD APPLICATION NODE (THE NODE FX,FY,FM ARE 
C APPLIED 

C DOF DEGREE OF FREEDOMS (UNKNOWN DISPLACEMENTS & SLOPES) 

C DVIBG. .. .DESIGN VARIABLE tl (BASE DIMENSION) INITIAL ESTIMATE 

C DVILO DESIGN VARIABLE #1 (BASE DIMENSION) LOV/ER SIDE CONSTRAINT 

C DVIUP DESIGN VARIABLE #1 (BASE DIMENSION) UPPER SIDE CONSTRAINT 

C DV2BG DESIGN VARIABLE #2 (SLOPE) INITIAL ESTIMATE 

C DV2L0 DESIGN VARIABLE t2 (SLOPE) LOV/ER CONSTRAINT 

C DV2UP DESIGN VARIABLE #2 (SLOPE) UPPER CONSTRAINT 

C EK 6X6 ELEMENT STIFFNESS MATRIX IN LOCAL X,Y COORDINATES 

C EKPR 6X6 ELEMENT STIFFNESS MATRIX IN ELEMENT LOCAL COORDINATES 

C ELEN LENGTH OF ELEMENT 

C F FORCE VECTOR OF SYSTEM 

C FA CONSTANT DISTRIBUTED LOAD OUTV/ARD FROM END TO END 

C FM CONCENTRATED MOMENT AT CLAN 

C FX CONCENTRATED LOAD IN X DIRECTION AT CLAN 

C FY CONCENTRATED LOAD IN Y DIRECTION AT CLAN 

C G THE ARRAY OF CONSTRAINT FUNCTIONS 

C GAMMA 6X6 ELEMENT TRANSFORMATION MATRIX 

C CK (NDOF)X(NDOF) GLOBAL STIFFNESS MATRIX 

C HGT CONSTANT DEPTH OF CROSS-SECTION 

C INFO ADS PARAMETER USED TO SIGNAL THAT THE OPT IS COMPLETE 

C IPRINT...ADS PARAMETER USED SELECT THE DATA OUTPUT FORMAT 

C ITERATE.. THE NUMBER OF TIMES ADS IS TO EE RELOADED V/ITH THE 

C PRECEEDING DATA 

C IV/K ADS INTERNAL WORK SPACE ARRAY 

C NCON NUMBER OF DESIGN CONSTRAINTS 

C NDOF NUMBER OF DEGREES OF FREEDOM 

C NDV NUMBER OF DESIGN VARIABLES 

C NEL NUMBER OF ELEMENTS 

C NRIWK....ADS INTERNAL WORK SPACE ARRAY DIMENSION 

C NRV/K ADS INTERNAL WORK SPACE ARRAY DIMENSION 

C NSNP NUMBER OF SYSTEM NODAL POINTS 

C OBJ THE OBJECTIVE FUNCTION OF THE OPTIMIZATION 

C OPTDCS... OPTIMIZATION DECISION TO OPTIMIZE THE PROBLEM OR NOT 
C PI.. .P5.. PARAMETER DIMENSION CORRESPONDING TO THE NEL, NSNP, NCON, 

C NDOF, AND NDV RESPECTIVELY 

C PRCSN....THE PRECISION DESIRED TO SOLVE THE FEM SYSTEM OF EQUATIONS 
C SIGMA_B..THE ELEMENTAL NORMAL STRESS DUE TO BENDING 

C SIGMA_N..THE ELEMENTAL NORMAL STRESS DUE TO AXIAL FORCES 

C SIGMA_T..THE MAXIMUM TOTAL STRESS IN EACH ELEMENT 

C SX GLOBAL HORIZONTAL COORDINATE 

C SY GLOBAL VERTICAL COORDINATE 

C U THE "DISPLACEMENT” VECTOR OF THE SYSTEM OF LINEAR EQUATIONS 

C VLB ADS ARRAY CONTAINING UPPER SIDE CONSTRAINTS 

C VUB ADS ARRAY CONTAINING LOWER SIDE CONSTRAINTS 

C WK ADS INTERNAL WORK AREA 

C X ADS ARRAY CONTAINING THE VALUES OF THE DESIGN VARIABLES 

C YIELD YIELD STRENGTH OF THE ARCH MATERIAL 

C YOUNG. .. .YOUNG'S MODULUS OF THE ARCH MATERIAL 

Cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C23A567 DECLARE THE VARIABLES 

INCLUDE 'ARCHCOM FORTRAN* 

C 



75 



oo oo oo oooooo oo o oo ooo oo 



FILE: ARCHOPT FORTRAN A 



C 

C 



....read the input parameters 

0PEN(8, FILE=*ARCHIN*, STATUS= *OLD » ) 



READ(8,)( ) 
READ(8,x) 
READ(8,x ) 
READ(8,x ) 
READ(8,x) 
READ(8,x) 
READ(8,X) 
READ(8,x ) 
READ(8,x) 
READ(8,x) 



L,H,HGT 
YOUNG, YIELD 

NEL,ISTRAT,IOPT,IONED,IPRINT,IGRAD 

DV1BG,DV1L0,DV1UP 

DV2BG,DV2L0,DV2UP 

CLAN,FX,FY 

FM,FA,OPTDCS 

ITERATE, PRCSN,BX1,BY1,BM1 

BX2,BY2,BM2 

LABEL 



....define constants 
NSNP = NEL + 1 
NDOF = 3XNSNP 
NOON = NSNP 

NDV = NSNP + (NEL - 1) 



OPTIMIZE THE PROBLEM 

CALL 0PTIMI2ATI0N_T00L 

COMPILE AND FORMAT THE OUTPUT 

CALL ARCH^OUTPUT 



END 



SUBROUTINE OPTIMIZATION^TOOL 



THIS SUBROUTINE DIRECTS THE PROGRAM FLOW OPTIMIZATION DECISION 
I.E., OPTIMIZE THE PROBLEM OR NOT. IT ALSO SERVES TO SET UP & 
EXECUTE THE ADS OPTIMIZATION SOFTWARE. 



....declare the variables 

INCLUDE *ARCHCOM FORTRAN' 

INTEGER I 
DO 100 1=1, NSNP 

EASE(I) = DVIBG 
BASEL(I) = DVILO 
BASEUm = DVIUP 
100 CONTINUE 

DO 150 I=1,NEL-1 

SLPm = DV2BG 
SLPL(I) = DV2L0 
SLPU(I) = DV2UP 
150 CONTINUE 

SLP(NEL)=DV2BG 

SLPL(NEL)=DV2L0 

SLPU(NEL)=DV2UP 

COMBINE BASE AND SLP ARRAYS INTO DESIGN ARRAY. 

DO 200 1=1, NSNP 

X(I) = BASE(I) 

VLB(I) = BASEL(I) 

VUB(I) = BASEU(I) 

200 CONTINUE 

DO 250 J=NSNP+1,NDV 

X(J) = SLP(J-NSNP) 

VLB( J)=SLPL( J-NSNP) 

VUB( J)=SLPU( J-NSNP) 

250 CONTINUE 

MAKE OPTIMIZATION DECISION 

IF (OPTDCS .NE. 1) THEN 
CALL EVAL 
RETURN 
ENDIF 



DEFINE THE SIZE OF THE WORK ARRAYS FOR ADS 
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NRA=65 
NC0LA=97 
NRWK = AOOOO 
NRIWK = 2000 

C ....ready to optimize... 

DO 280 I=1,NSNP 
1DG(I)=0 
280 CONTINUE 
INFO =-2 

CALL ADS ( INFO, ISTRAT, lOPT, lONED, IPRINT, IGRAD,NDV,NCON, X, VLB, VUB, 
: OBJ, G, IDG, NGT,IC,DF, A, NRA,NCOLA,WK, NRWK, IV/K, NRIWK) 

IV/K(2)=0 



300 CALL ADS ( INFO, ISTRAT, lOPT, lONED, IPRINT, IGRAD, NDV,NCON, X, VLB, VUB, 
: OBJ, G, IDG, NGT,IC,DF, A, NRA,NCOLA,WK, NRWK, IWK, NRIWK) 

C 

C ....evaluate the objective function and constraints 

IF (INFO .NE. 0) THEN 
CALL EVAL 
GOTO 300 

ENDIF 

END 

CMXKKMKXXXXKMKXXKXMMiOOOOCXXMKKKXXKlOnCXKXXMXinnnCKXKMMKMXXKMMMMXMXKXXMMX 

C 

SUBROUTINE EVAL 

C ============================z==========r===z======================== 

C THIS SUBROUTINE IS USED TO EVALUATE THE OBJECTIVE FUNCTION, 

C CONSTRAINT FUNCTIONS, AND SIDE CONSTRAINTS OF THE OPTIMIZATION 
C PROBLEM. 

C ============z===========z===z==z=zz========z===========r=========r== 



C ....declare the variables 

INCLUDE 'ARCHCOM FORTRAN’ 

INTEGER I,J,K 
REAL PI 

PARAMETER ( PI=3. 1A1593 ) 

C 

C SEPARATE THE DESIGN ARRAY INTO EASE AND SLP ARRAYS 

DO 50 I=1,NSNP 

BASE(I) = X(I) 

50 CONTINUE 

DO 75 I=1,NDV-NSNP 

SLP(I) = XCI+NSNP) 

75 CONTINUE 
DX=L/NEL 
SY(1)=0.0 
DO 80 I=1,NEL-1 

SY(I + l)=DX3tSLP(I) + SY(I) 

80 CONTINUE 

SLP(NEL)=(H - SY(NEL))/DX 
C 

C ....calculate the objective function 

OBJ = 0.0 

c 

DO 100 1=1, NEL 

BAVE(I) = (BASE(I)+BASE(I+l))/2.0 
OBJ = OBJ + BAVE(I)XSQRT(1 + SLP(I)*«2) 

100 CONTINUE 

C CALCULATE ALPHA( I ) ,BETA( I ) , AND ELEN(I) 

DO 150 1=1, NEL 

ALPHA( I ) =ATAN2 ( SLP ( I )xDX , DX ) 

BETA(I)=(PI/2) - ALPHA(I) 

ELEN(I)=SQRT(DXxx2 + (SLP( I )XDX )xx2 ) 

150 CONTINUE 
C 

C ....determine the design constraints.... 

CALL ARCH^STRESS 
C 

DO 230 I=1,NSNP 

G(I) = (SIGMA_T( D/YIELD - 1.0) 

230 CONTINUE 
END 

CKKKXXKKXXXKMXMMMXKKKKMXMMXXXKXMXKXXKKXXKKKKKMMMKKKKXKXKMMXXXKMMXXKKXK 



77 



oooooo 



file: archopt fortran a 



SUBROUTINE ARCH.STRESS 



THIS SUBROUTINE IS USED TO PERFORM THE FINITE ELEMENT ANALYSIS 
OF THE STRESSES DEVELOPED IN AN ARCH OR BEAM FOR A GIVEN LOAD- 
ING. 



....declare the variables 

INCLUDE ’ARCHCOM FORTRAN' 

INTEGER IPVT(99) 

REAL F(P4J 

REALM 8 BK ( PA , PA ) , BF ( PA ) , BU ( P4 ) , FAC ( 980 1 ) » WORK ( 99 ) 



C 

C ....form the element and system matrices 

CALL FORM 
C 

C FORM THE FORCE VECTOR, F 

CALL FORCE.VECTOR (NEL,NDOF, ELEN, ALPHA, BETA, FA, F) 

C 

C ....SET THE BOUNDARY CONDITIONS AND LOADS 

CALL ENDARY ( NDOF, 6K, CLAN, FX, FY, FM, F, BXl, BY1,BM1,BX2,BY2, BM2 ) 

C 

C ....SOLVE THE SYSTEM OF EQUATIONS 

IF (PRCSN .EQ. 2) THEN 

C CHANGE GK AND F ARRAYS TO DOUBLE PRECISION 

CALL UPSCALE (NDOF,GK, F,EK,BF) 

C SOLVE THE SYSTEM OF EQUATIONS 

CALL DL2ARG ( NDOF,BK, PA, BF, 1 , EU, FAC , IPVT, V/ORK ) 

C CHANGE EU ARRAY TO SINGLE PRECISION 

CALL DOV/NSCALE (NDOF,EU,U) 

ELSE 

C SOLVE THE SYSTEM OF EQUATIONS 

CALL L2ARG ( NDOF, GK, PA, F, 1, U, FAC, IPVT, V^ORK ) 

ENDIF 

C 

C ....determine the stress distribution 

CALL STRESS 
C 

C RETURN 
END 

Cmmmmmmmmxmxxxxxxmmmxxmmmmmmmmmmxxxmmmxxmxmxxxxxxxmxmmmmmxmmmmmmmmmmmm 
SUBROUTINE FORM 

C r=====r===rr==s==============================r=====z=====r========== 

C This subroutine is used to construct the global stiffness mat- 

C RIX FOR THE ARCH PROBLEM. 

C ==================================================================== 

C ....declare the variables 

INCLUDE 'ARCHCOM FORTRAN' 

INTEGER lEL, I, J,K,II, JJ,KK,III, JJJ 

REAL C1,C2,C3,CA,C5,CA,CB,EK(P1,6,61,GAMMA(6,6),EKGA(6,6), 

: GAEKGA(6,6),BH,BH3 

C 

DO 120 IEL::1,NEL 

C ....define the constants Cx 

Cl = YOUNG/ELENdEL) 

C2 = (1.0/ELEN(IEL))m*2.0 
C3 = (1.0)/(2.0kELEN(IEU) 

CA = ll.OJ/3.0 
C5 = (l.OJ/6.0 
C 

C ....initialize the v/ork arrays 

DO 100 I = 1,6 
DO 90 J= 1,6 

EKPR(IEL,I,J) = 0.0 
GAMMA(I,J) = 0.0 
EKGA(I,J) = 0.0 
GAEKGA(I,J) = 0.0 
EK(IEL,I,J) = 0.0 
90 CONTINUE 

100 CONTINUE 

C 

C ....calculate the area and inertia terms 
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C 

C 



C 

C 



BH = BAVE(IEL)«HGT 

EH3 = BAVE(ZEL)X(HGTXX3.0) 



.... determine 
EKPR(IEL,1,1) 
EKPR(lEL,lrA) 
EKPR(IEL,2^2) 
EKPR(IEL,2,3) 
EKPR(IEL,2,5) 
EKPRC IEL,2,6) 
EKPR(IEL,3,2) 
EKPR( IEL,3f3) 
EKPR(IEL,3,5) 
EKPR(1EL»3»6) 
EKPR(IEL/A,1) 
EKPR(IEL,A,A) 
EKPR(IEL,5,2) 
EKPR(IEL,5,3) 
EKPR(IEL,5,5) 
EKPR(IEL,5,6) 
EKPR(IEL,6,2) 
EKPR(IEL,6,3) 
EKPR(IEL,6,5) 
EKPR(IEL,6,6) 



the EKPR matrix, 
= ClxBH 
= -ClxBH 
= C1«C2»EH3 

= C1XC3XBH3 
= -C1«C2*BH3 
= C1«C3*BH3 

= C1*C3«EH3 

= C1XCAXBH3 
= -C1*C3*BH3 
= C1mC 5«BH3 
= -ClxBH 
= C1«BH 
= -C1XC2XBH3 
= -C1 mC 33(BH3 
= C1XC2XBH3 
= -C1XC3XBH3 
= C1XC3XEH3 
= ClxC5*EH3 
= -C1XC3XBH3 
= ClxCA3<BH3 



....determine the GAMMA matrix 
CA = COS(ALPHAdEL)) 

CB = COS(BETA(IEU) 
GAMMA(1,1) = CA 
GAMMA(1.2) = CE 
GAMMA(2,1) = -CB 
GAMMA(2,2) = CA 
GAMMA(3,3) = 1.0 

GAMMA(A,A) = CA 
GAMMA(A,5) = CB 
CAMMA(5,A) = -CB 
GAMMA(5,5) = CA 
GAMMA(6,6) = 1.0 



C ....determine the EKGA array 

DO 220 I = 1,6 
DO 215 J = 1,6 
BO 210 K = 1,6 

EKGA(I,J) = EKCA(I,J) + EKPR( lEL, I, K)»GAMMA( K, J ) 
210 CONTINUE 

215 CONTINUE 

220 CONTINUE 

C 

C ....determine the GAEKGA array 

DO 2A0 I = 1,6 
DO 235 J = 1,6 
DO 230 K s 1,6 

GAEKGA(I,J) = GAEKGA(I,J)+GAMMA(K,I)»EKGA(K,J) 
230 CONTINUE 

235 CONTINUE 

2A0 CONTINUE 

C 

C ....copy the GAEKGA array into the EK array 

DO 260 I = 1,6 
DO 250 J = 1,6 

EK(IEL,I,J) = GAEKGA(I,J) 

250 CONTINUE 

260 CONTINUE 

120 CONTINUE 
C 

C ....initialize the GK array 

DO 150 1=1, NDOF 
DO lAO J = 1, NDOF 
GK(I,J) = 0.0 
lAO CONTINUE 

150 CONTINUE 



79 



FILE: ARCHOPT FORTRAN A 



C 

C ....construct the GK matrix 

DO 300 lEL = 1, NEL 

II = 3x(IEL-l) 

DO 290 J = 1, 6 

JJ = II + J 
DO 280 K = 1, 6 

KK = II + K 

GKCJJrKK) = GKC JJ,KK)+EKCIEL, J,K) 

280 CONTINUE 

290 CONTINUE 

300 CONTINUE 
C 

C RETURN 

END 

SUBROUTINE FORCE_VECTOR ( NEL, NDOF, ELEN , ALPHA, BETA, FA, F ) 

C ==================================================================== 

C This subroutine is used to construct the force vector for the 

C FEM PROBLEM SPECIFIED. 

C ==================================================================== 

C ....DECLARE THE VARIABLES 

INTEGER NEL, NDOF, I, II, 12, 13, PI, PA 
C 

PARAMETER( Pl=32, PA=99 ) 

C 

REAL ELEN(P1),ALPHA(P1),BETA(P1),FA,F(PA) 

C 

C ...FORM THE F-VECTOR 

FCll = (ELEN(l]/2.0 J x ( -COS (BETAf 1 ) )) 

F(2) = (ELEN(l)/2.0) x { COS ( ALPHAt 1 ) ) ) 

FC3) = 0.0 
C 

DO 100 1=2, NEL 

11 = (I-l)x3 + 1 

12 = (i-l)x3 + 2 

13 = (i-l)x3 + 3 
C 

F(I1) = (ELEN(I)/2.0)X(-COS(BETACI) ) ) 

: +tELEN(I-l)/2.0)xC-COS(EETA{I-l) ) ) 

F(I2) = tELENfI)/2.0)xlCOS(ALPHACI) n 
: + 1 ELEN ( I-l )/2. OJxt COS (ALPHA (I-l) ) ) 

p(I3) = 0.0 
100 CONTINUE 
C 

F(NDOF-2) = (ELEN(NEL)/2.0)x(-COS(BETA(NEU ) ) 

F(NDOF-l) = (ELEN(NEL)/2.0)x(COS(ALPHA(NEU ) ) 

F(NDOF) =0.0 
C 

C SCALE THE F-VECTOR BY FA 

DO 200 1=1, NDOF 
F(I) = FAXF(I) 

200 CONTINUE 
C 

C RETURN 

END 

Cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SUBROUTINE BNDARY (NDOF, GK, CLAN, FX, FY, FM, F,BX1,BY1,BM1,BX2, 



C This subroutine is used to impose the boundary conditions upon 

C THE GLOBAL STIFFNESS MATRIX AND FORCE VECTOR. 

C ==================================================================== 

C ....declare the variables 

INTEGER ND0F,BX1,BY1,BM1,BX2,BY2,BM2,CLAN,I,N,I1,I2,I3,PA 
PARAMETER(PA=99) 

REAL GK(PA,PA),FX,FY,FM,F(PA) 

C 

C ....invoke the essential boundary conditions 

IF (BXl .EQ. 1) THEN 
CALL IMPOSEBC (NDOF, GK, 1, F ) 

ENDIF 
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C 

IF (BYl .EQ. 1) THEN 

CALL IMPOSEBC ( NDOF, GK, 2, F ) 

ENDIF 

C 

IF (EMI .EQ. 1) THEN 

CALL IMPOSEBC ( NDOF, GK,3, F ) 

ENDIF 

C 

IF (BX2 .EQ. 1) THEN 
N=ND0F-2 

CALL IMPOSEBC ( NDOF, GK, N, F ) 

ENDIF 

C 

IF (BY2 .EQ. 1) THEN 
N=NDOF-l 

CALL IMPOSEBC ( NDOF, GK, N, F ) 

ENDIF 

C 

IF (BM2 .EQ. 1) THEN 

CALL IMPOSEBC ( NDOF, GK, NDOF, F ) 

ENDIF 

C 

C ADD THE CONCENTRATED LOAD TO THE FORCE VECTOR 

Il=(CLAN-l)x3+l 

I2=(CLAN-l)«3+2 

I3=(CLAN-l)»3+3 

C 

F(I1)=F(I1)+FX 

F(I2)=F(I2)+FY 

F(I3)=F(I3)+FM 

C 

C RETURN 

END 

CXXKXXXXXKXXXKXXMXXXXXXMXXXMXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXMMXMXXMXXX 
SUBROUTINE IMPOSEBC (NDOF,GK,N,F) 

C ==================================================================== 

C This subroutine is used to do the redundant leg work of impos- 

C ING THE BOUNDARY CONDITIONS. 

C ==================================================================== 

C DECLARE THE VARIABLES 

INTEGER NDOF,N,I,PA 
PARAMETER(P4=99) 

REAL GK(P^,P4),F(PA) 

C 

C ....IMPOSE THE BOUNDARY CONDITION ON THE GK AND F ARRAYS 

DO 100 I=l,NDOF 
GK(N,I) = 0.0 
100 CONTINUE 

GK(N,N) = 1.0 
F(N) = 0.0 
C 

C RETURN 

END 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
SUBROUTINE UPSCALE (NDOF, GK, F,EK,BF) 

C ==================================================================== 

C This subroutine is used to change the stiffness matrix & force 

C VECTOR FROM SINGLE PRECISION TO DOUBLE PRECISION IN ORDER TO 

C SOLVE THE LINEAR SYSTEM OF EQUATIONS IN DOUBLE PRECISION. 

C ==================================================================== 

C ....DECLARE THE VARIABLES 

INTEGER NDOF,I,J,P^» 

PARAMETER (PA=99) 

REAL GK(PA,P4),F(P^i) 

REALxS BK(P^,PA),BF(PA) 

....GENERATE THE DOUBLE PRECISION COMPLIMENTS OF GK AND F 

DO 110 I=l,NDOF 
DO 100 J=l,NDOF 
BK(I,J) = GK(I,J) 

100 CONTINUE 
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BF(I) = F(I) 

110 CONTINUE 
C 

C RETURN 

END 

CXXXXXXXKXXXXXKXMXXKXXXMXXXXMXXXMXXXXXXKKKMKKKXKXXXXXKKKXXXXXXKXXXXXXXK 

SUBROUTINE DOWNSCALEt NDOF,BU, U ) 

C =====rr==================r==r=r=z================rr=====zz========== 

C This subroutine is used to do down scale the double precision 

C SOLUTION OF THE LINEAR SYSTEM OF EQUATIONS BACK TO SINGLE PRE- 

C CISION, ADS COULD HAVE PROBLEMS WITH DOUBLE PRECISION NUMBERS! 

C =================================================== ============r=zzz 

C ....declare the variables 

INTEGER NDOF,I,PA 
PARAMETER ( PA=99 ) 

REAL U(PA) 

REALX8 BU(PA) 

C 

C GENERATE THE SINGLE PRECISION COMPLIMENT OF BU 

DO 100 I=l,NDOF 

U(i) = BU(i) 

100 CONTINUE 
C 

C RETURN 
END 

Cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SUBROUTINE STRESS 

C =================================================================== 

C THIS SUBROUTINE COMPUTES THE STRESS AT EACH NODAL POINT. 

C =================================================================== 

C ....declarations 

INCLUDE 'ARCHCOM FORTRAN' 

INTEGER II, 12, 13, lA, 15, 16, I, Jr K 

REAL CA1,CB1,K1,K2,FPR(PA,6),UPR(6),N0RM1,N0RM2, 

BEND1,BEND2 



C 



c 



c 



250 

C 



350 

300 

100 



....determine local 
DO 100 1=1, NEL 

I1=(I-1)X3+1 

I2=(i-l)x3+2 

I3=(i-l)x3+3 

IA=(i)x3+l 

I5=(i)x3+2 

I6=( i)x3+3 



forces from stiffness 



and displacement.... 



C31= COS(EETA(D) 
CA1= COS(ALPHAd)) 



UPR(1)= UdDxCAl + 
UPR(2)= -UdDxCBl + 
UPR(3)= Ud3) 
UPR(A)= UdAJxCAl + 
UPR(5)= -UdA)xCBl + 
UPR(6)= Ud6) 



Ud2)xCBl 

Ud2)xCAl 

Ud5)xCBl 

Ud5)xCAl 



DO 250 JJ=1,6 

FPRd,JJ)= 0.0 
CONTINUE 



DO 300 J=l,6 
DO 350 K=l,6 

FPRd,J)= FPRd,J) + EKPRd, J,K)XUPR(K) 
CONTINUE 
CONTINUE 
CONTINUE 

....determine the bending and normal stresses 

SIGMA.Nd) = ABS(FPRd,l)xd.O/(BASEd)xHGT))) 
SIGMA_Bd) = ABS(FPRd,3)x(6.0/(BASE(l)x(HGTXx2.0)))) 
SIGMA_Td) = SIGMA_Bd) + SIGMA_Nd) 

DO AOO i=2,NEL 

K1 = 1.0/(BASEd)xHGT) 

K2 = 6.0/(BASEd)x(HGTxx2.0)) 
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NORMl = ABS(FPR(i,l)«Kl) 
N0RM2 = ABS(FPR(i-l,4)«Kl) 

BENDl = ABS(FPR(i,3)«K2) 
BEND2 = ABS(FPR(i-l,6)«K2) 



SIGMA_N(i) = (N0RMl+N0RM2)/2.0 
SIGMA_B(i) = (BENDl+BEND2)/2.0 



SIGMA_T(i) 



SIGMA_B(i) + SIGMA_N(i) 



AOO CONTINUE 



SIGMA_N(NSNP) = ABS( FPR( NEL, 4 )« ( 1. 0/(BASE (NSNP )«HGT ) ) ) 
SIGMA_B(NSNP) = ABS( FPR(NEL,6 )« 

: (6.0/(BASE(NSNP)«(HGT«m 2.0) ) n 

SIGMA^T(NSNP) = SIGMA.B (NSNP 3 + SIGMA_N (NSNP ) 

C 

C RETURN 

END 



Cmmmxmxxxxxxxxmxxxxmxmxxmxxmxxxxxmmmmxmxmmmmmmxmxxxmxmmmmxmxmxxmmmmmmm 
SUBROUTINE ARCH.OUTPUT 



THIS SUBROUTINE FORMATS THE FINAL RESULTS AND OUTPUT OF THE 
OPTIMIZATION PROBLEM AND STORES IT IN A FILE NAtlED ARCH_OUT 



DECLARE VARIABLES 

INCLUDE 'ARCHCOM FORTRAN' 
REAL VOL, VOLUME 



OPEN OUTPUT FILE AND V/RITE HEADER 

0PENC9, FILE='ARCHOUT', STATUS= ' OLD ' ) 

C 

WRITE(9,100) LABEL 

V/RITE(9,100) ' OPTIMIZATION SOLUTION' 

V^RITE(9, 105) ' — 

. I 

100 * FORMAT (/5X, A) 

105 FORMAT (5X, A) 

C 

C SECTION "A" 

V/RITE(9,100) ' A) PROBLEM PARAMETERS:' 

URITE(9,110) ' HORIZ. SPAN:', L, ' YOUNGS MODULUS: ', YOUNG 
WRITE(9,110) ' VERT. SPAN :', H, ' YIELD STRENGTH: », YIELD 
WRITE(9, 115) ' NO OF DESIGN VAR:', NDV, ' NO OF ELEMENTS: ' ,NEL 
110 F0RMAT(8X,A,F12.3,T38,A,F12. 1) 

115 FORMAT(8X,A,I7,T38,A,I10) 

C 

C SECTION "B" 

WRITE(9,100) ' B) DERIVED CONSTANTS:' 

WRITE(9,120) ' NO OF SYSTEM NODAL POINTS. ,NSNP 

WRITE(9,120) ' NO OF DEGREES OF FREEDOM ',NDOF 

V/RITE(9,125) ' HORIZ. LENGTH PER ELEMENT. DX 
120 F0RMAT(8X,A,I6) 

125 F0RMAT(8X,A,F12.4) 

C 

C ....section "C" 



V/RITE(9,100) 
V/RITE(9,125) 
WRITE (9, 125) 


' C) STRUCTURE LOADING:' 

» FX 


* »FX 


» FY 


* f FY 


V/RITE(9,125) 

WRITE(9,125) 

WRITE(9,120) 


'FM 

'FA 

'CONCENTRATED LOAD AT NODE... 


',FM 

SFA 

',CLAN 



C . . . .section "D" 

SX(1)=0.0 
SY(1)=0.0 
DO 150 1=1, NEL 
SX(I+l)=lxDX 

SY(I+1)=DXXSLP(I) + SY(I) 
150 CONTINUE 
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WRITE(9,100) * D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION:* 
WRITE(9,210) *NODE*, *X-COORD*, 'Y-COORD*, 'BASE*, 'AREA' 

C 

210 F0RMAT(8X,A,T21,A,T36,A,T49,A,T62,A) 

220 FORMAT(8X,IA,T17,F10.5,T32,F10.S,TA8,F8.5,T60,F8.5) 

VOLUME =0.0 
C 

DO 300 I=1,NSNP 

AREA = BASE(I)xHGT 

WRITE (9, 220) I , SX ( I ) , SY( I ) ,BASE( I ) , AREA 
300 CONTINUE 
C 

C . . . .section "E" 

VOL = HGTXDXXOBJ 

WRITE(9,100) ' E) OBJECTIVE FUNCTION:* 

C 

WRITE(9,310) * TOTAL STRUCTURE VOLUME:*, VOL 
310 F0RMAT(/12X,A,F12.6/) 

C 

WRITE(9,330) 'NODE* , 'NORMAL STRESS' , 'BENDING STRESS* , 'TOTAL* 

DO 320 I=1,NSNP 

WRITE(9,3A0) i,SIGMA_N( i ) ,SIGMA_B( i ) , SIGMA_T( i ) 

320 CONTINUE 

330 F0RMAT(8X,A,T18,A,T35,A,T57,A) 

3A0 F0RMAT(8X,IA,T15,F1A. 1,T32,F1A.1,T49,F1A. 1) 

C 

C . . . .section "F" 

V/RITE(9, 100) * F) BOUNDARY CONDITIONS:' 

WRITE(9,A10) 'NODE*, 'X-DISPL' , 'Y-DISPL* , 'SLOPE' 

WRITE(9,A30) 1 , EXl, BYl , EMI 
WRITE(9,A30) NEL+1,BX2,BY2,BM2 
C 

C SECTION "G" 

WRITE(9,100) * G) SOLUTION VECTOR:' 

WRITE(9,A10) 'NODE', 'X-DISPL*, 'Y-DISPL*, 'SLOPE* 

DO AOO I=1,NSNP 
I1=(I-1)3(3+1 
I2=(I-l)x3+2 
I3=(I-l)*3+3 

WRITE(9,A20) I , U( I 1 ) ,U( 12 ) , U( 13 ) 

AOO CONTINUE 

A 10 FORMATCT9,A,T17,A,T31,A,TA6,A) 

A20 F0RMATC7X,I5,3E1A.6) 

A30 FORMAT(7X,I5,T20,IA,T3A,IA,TA8,IA) 

C 

C RETURN 

END 
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FILE: FILE 



CASEl 



A appendix C: CASE STUDIES 



OPTIMIZATION NO.l 
OPTIMIZATION SOLUTION 



A I PROBLEM PARAMETERS: 

HORIZ. SPAN: 32.000 

VERT. SPAN : 32.000 

NO OF DESIGN VAR: 2A 

B) DERIVED CONSTANTS: 

NO OF SYSTEM NODAL POINTS... 

NO OF DEGREES OF FREEDOM 

HORIZ. LENGTH PER ELEMENT... 

C) STRUCTURE LOADING: 

FX 

FY 

FM 

FA 

CONCENTRATED LOAD AT NODE... 



YOUNGS MODULUS: 30000000.0 
YIELD STRENGTH: 52000.0 
NO OF ELEMENTS: 12 



13 

39 

2.6667 



0.0000 

0.0000 

0.0000 

- 100.0000 

7 



D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 



NODE 


X-COORD 


Y-COORD 


BASE 


AREA 


1 


0.00000 


0.00000 


0.15000 


0.22500 


2 


2.66667 


A. 30175 


0. 15000 


0.22500 


3 


5.33333 


8.23110 


0.1508A 


0.22627 


A 


8.00000 


11.7A150 


0.15022 


0.22533 


5 


10.66666 


lA. 87895 


0.15092 


0.22633 


6 


13.33333 


17.7212A 


0. 1501A 


0.22522 


7 


16.00000 


20.31323 


0.15000 


0.22500 


8 


18.66666 


22.73A09 


0. 15000 


0.22500 


9 


21.33333 


2A. 91997 


0. 15000 


0.22500 


10 


23.99998 


26.86A90 


0.15000 


0.22500 


11 


26.66666 


28.66069 


0.15000 


0.22500 


12 


29.33331 


30.59526 


0. 15000 


0.22500 


13 


31.99998 


31.99998 


0. 15000 


0.22500 



E) OBJECTIVE function: 



TOTAL STRUCTURE VOLUME: 10.32A365 



NODE 


NORMAL STRESS 


BENDING STRESS 


TOTAL 


1 


36039.3 


0.1 


36039.5 


2 


36031.7 


15159.3 


51191.0 


3 


35823.3 


16175.0 


51998.3 


A 


35975.3 


15313.8 


51289.2 


5 


35809.8 


15699.2 


51509. 1 


6 


3599A.7 


16009.5 


5200A.3 


7 


36033.5 


15762. A 


51796.0 


8 


360AA.7 


9766.7 


A5811.A 


9 


360A9.8 


7516.0 


A3565.8 


10 


360A1. A 


12A57.1 


A8A98.6 


11 


36016.8 


16026.7 


520A3.5 


12 


36039.7 


15969.3 


52009.0 


13 


36081. A 


0.0 


36081. A 



F) BOUNDARY CONDITIONS: 
NODE X-DISPL 


Y-DISPL 


1 


1 


1 


13 


1 


1 



SLOPE 

0 

0 



G) SOLUTION VECTOR! 
NODE X-DISPL 

1 O.OOOOOOE+00 

2 0.509330E-01 

3 0.S63026E-01 
A 0.106222E+00 

5 0.11AA65E+00 

6 0.113685E+00 

7 0.105703E+00 



Y-DISPL 

O.OOOOOOE'i'OO 

-0.387271E-01 

-0.696029E-01 

-0.913596E-01 

-0.10A831E+00 

-0.110A95E+00 

-0.108688E+00 



SLOPE 

-0.131531E-01 

-0.11AA81E-01 

-0.81A113E-02 

-0.505612E-02 

-0.221819E-02 

0.528031E-03 

0.315360E-02 
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A 



8 0.921987E-01 -0. 1002^i9E+00 

9 0,7612A2E-01 -0. 871751E-01 

10 0.589296E-01 -0. 703301E-01 

11 0.3963A1E-01 -0. A85910E-01 

12 0.16288AE-01 -0. 231^»30E-01 

13 O.OOOOOOE+00 O.OOOOOOE+00 



0.519679E-02 
0.652113E-02 
0.798611E-02 
0. 100210E-01 
0. 100252E-01 
0.895565E-02 
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CASEIA A 



OPTIMIZATION NO.XA 
OPTIMIZATION SOLUTION 



A) PROBLEM parameters: 

HORIZ. SPAN: 32.000 

VERT. SPAN : 32.000 

NO OF DESIGN VAR: 2A 

E) DERIVED CONSTANTS: 

NO OF SYSTEM NODAL POINTS.. 
NO OF DECREES OF FREEDOM... 
HORIZ. LENGTH PER ELEMENT.. 

C) STRUCTURE LOADING: 

FX 

FY 

FM 

FA 

CONCENTRATED LOAD AT NODE... 



YOUNGS MODULUS: 30000000.0 
YIELD STRENGTH: 52000.0 
NO OF ELEMENTS: 12 



13 

39 

2.6667 



0.0000 

0.0000 

0.0000 

- 100.0000 

7 



D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 



NODE 


X-COORD 


Y-COORD 


BASE 


AREA 


1 


0.00000 


0.00000 


0.79198 


1. 18797 


2 


2.66667 


2.66933 


0.78A08 


1. 17612 


3 


5.33333 


5.33867 


0.78335 


1.17503 


A 


8.00000 


8.00800 


0.986AA 


1.A7966 


5 


10.66666 


10.67200 


1.16866 


1.75299 


6 


13.33333 


13.33852 


1.27656 


1.91A8<» 


7 


16.00000 


16.00357 


1.30861 


1.96291 


8 


18.66666 


18.66756 


1.27A57 


1.91186 


9 


21.33333 


21.3315A 


1.16A62 


1.7A693 


10 


23.99998 


23.99666 


0.98731 


1.A8096 


11 


26.66666 


26.66280 


0.7832A 


1. 17A86 


12 


29.33331 


29.32883 


0.78168 


1.17252 


13 


31.99998 


31.99998 


0.78982 


1.18A73 



E) OBJECTIVE FUNCTION: 



TOTAL STRUCTURE VOLUME: 68.3867A9 



NODE 


NORMAL STRESS 


BENDING STRESS 


TOTAL 


1 


91.5 


0.0 


91.5 


2 


92.6 


26676. A 


26769.0 


3 


92.2 


AS556.6 


A86A8.8 


A 


72.9 


52068.1 


521A0.9 


5 


61.6 


5209A.5 


52156. 1 


6 


56.0 


52183.3 


52239. A 


7 


5A.3 


52353.9 


52A08.3 


8 


56.5 


52263.8 


52320.3 


9 


62.2 


52273.6 


52335.8 


10 


73.0 


5201A. 1 


52087. 1 


11 


92.3 


A856A.5 


A8656.7 


12 


91.5 


26770.9 


26862. A 


13 


89.5 


0.2 


89.8 


F) BOUNDARY 


CONDITIONS: 







NODE X-DISPL Y-DISPL SLOPE 

11 1 0 
13 1 1 0 



G) SOLUTION vector: 
NODE X-DISPL 

1 O.OOOOOOE-t-00 

2 0.113A23E+00 

3 0.215283E+00 
A 0.297166E+00 

5 0.355750E+00 

6 0.390952E+00 

7 O.A02699E+00 



Y-DISPL 
O.OOOOOOE+00 
-0. 113326E+00 
-0.215100E+00 
-0.296917E+00 
-0.355571E+00 
-0.390785E+00 
-O.A025A9E+00 



SLOPE 

-0.A32361E-01 

-0.A10106E-01 

-0.3A7033E-01 

-0.262325E-01 

-0.17507AE-01 

-0.876822E-02 

-0.976327E-05 
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8 

9 

10 

11 

12 

13 



0.390992E+00 

0.355829E+00 

0.297182E+00 

0.215351E+00 

0.113571E+00 

O.OOOOOOE+00 



-0.3908A0E+00 
-0.355653E+00 
-0.296983E+00 
-0.215151E+00 
-0. 11336AE+00 
O.OOOOOOE+00 



0.875361E-02 
0. 17S100E-01 
0.262^»91E-01 
0.3^7106E-01 
0.A10251E-01 
0.A32590E-01 
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A 



OPTIMIZATION NO. 2 
OPTIMIZATION SOLUTION 



A) PROBLEM PARAMETERS: 

HORIZ. SPAN: 32.000 

VERT. SPAN : 32.000 

NO OF DESIGN VAR: 2A 

B) DERIVED CONSTANTS: 

NO OF SYSTEM NODAL POINTS.. 
NO OF DECREES OF FREEDOM... 
HORIZ. LENGTH PER ELEMENT.. 

C) STRUCTURE LOADING: 

FX 

FY 

FM 

FA 

CONCENTRATED LOAD AT NODE... 



YOUNGS MODULUS: 30000000.0 
YIELD STRENGTH: 52000.0 
NO OF ELEMENTS: 12 



13 

39 

2.6667 



0.0000 

- 2000.0000 

0.0000 

0.0000 

7 



D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 



NODE 


X-COORD 


Y-COORD 


BASE 


AREA 


1 


0.00000 


0.00000 


0. 15000 


0.22500 


2 


2.66667 


3.00301 


0. 15000 


0.22500 


3 


5.33333 


6.30157 


0. 15000 


0.22500 


A 


8.00000 


9.60067 


0.15000 


0.22500 


5 


10.66666 


12.91733 


0.15000 


0.22500 


6 


13.33333 


16.27870 


0. 15000 


0.22500 


7 


16.00000 


19.31017 


0. 15000 


0.22500 


8 


18.66666 


21.33359 


0.15000 


0.22500 


9 


21.33333 


23.39781 


0. 15000 


0.22500 


10 


23.99998 


25. A6082 


0. 15000 


0.22500 


11 


26.66666 


27.5A510 


0.15000 


0.22500 


12 


29.33331 


29.7A9AA 


0. 15000 


0.22500 


13 


31.99998 


31.99998 


0. 15000 


0.22500 



E) OBJECTIVE FUNCTION: 



TOTAL STRUCTURE VOLUME: 10.239A57 



NODE 


NORMAL STRESS 


BENDING STRESS 


TOTAL 


1 


31082.2 


0.3 


31082.6 


2 


31095.7 


20958. 1 


52053.9 


3 


31109.7 


18671.2 


A9780.9 


A 


31109. A 


16336.5 


A7AA5.9 


5 


31107.6 


12623.3 


A3730.9 


6 


31097. A 


5386. A 


36A83.7 


7 


27975.7 


2A09A.0 


52069.8 


8 


2A863.7 


27275.3 


52139.0 


9 


2A863.7 


27252.7 


52116. A 


10 


2A863.8 


27320.7 


5218A.5 


11 


2A857.5 


25720.0 


50577.5 


12 


2A8A6.A 


1A678.1 


3952A.5 


13 


2A8A1.8 


0.0 


2A8A1.8 



F) BOUNDARY CONDITIONS: 
NODE X-DISPL 

1 1 

13 1 



Y-DISPL 

1 

1 



SLOPE 

0 

0 



G) SOLUTION VECTOR: 
NODE X-DISPL 

1 O.OOOOOOE+00 

2 0.A61328E-01 

3 0.866829E-01 
A 0.115632E+00 

5 0.13A663E+00 

6 0.1A6308E+00 

7 0.151127E+00 



Y-DISPL 
O.OOOOOOE+00 
-0.A65305E-01 
-0.8A9686E-01 
-0.11A025E+00 
-0. 13A988E+00 
-0.1A9906E+00 
-0.159717E+00 



SLOPE 

-0.169057E-01 

-0.150352E-01 

-0.112998E-01 

-0.799962E-02 

-0.526077E-02 

-0.35A366E-02 

-0.898687E-03 
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8 

9 

10 

11 

12 

13 



0.1A69A9E+00 -0. 158800E+00 
0.13AA88E+00 -0. 1A7268E+00 
0.113602E+00 -0.12A838E+00 
0.8A0616E-01 -0.915992E-01 
0.AA8603E-01 -0 . A867AAE-01 
O.OOOOOOE+00 O.OOOOOOE+00 



0.292259E-02 
0.700891E-02 
0.110977E-01 
0.150871E-01 
0. 181930E-01 
0. 193311E-01 
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CASE3 



A 



OPTIMIZATION NO. 3 
OPTIMIZATION SOLUTION 



A) PROBLEM PARAMETERS: 

HORIZ. SPAN: 32.000 

VERT. SPAN : 32.000 

NO OF DESIGN VAR: 2A 

B) DERIVED CONSTANTS: 

NO OF SYSTEM NODAL POINTS.. 
NO OF DEGREES OF FREEDOM... 
HORIZ. LENGTH PER ELEMENT.. 

C) STRUCTURE LOADING: 

FX 

FY 

FM 

FA 

CONCENTRATED LOAD AT NODE... 



YOUNGS MODULUS: 30000000.0 
YIELD strength: 52000.0 
NO OF elements: 12 



13 

39 

2.6667 



2000.0000 

- 2000.0000 

0.0000 

0.0000 

7 



D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 



NODE 


X-COORD 


Y-COORD 


BASE 


AREA 


1 


0.00000 


0.00000 


0.15030 


0.22545 


2 


2.66667 


3.57692 


0.15037 


0.22556 


3 


5.33333 


7. 11881 


0.15045 


0.22568 


4 


8.00000 


11. 13781 


0.15079 


0.22618 


5 


10.66666 


14.64122 


0.15017 


0.22526 


6 


13.33333 


18.30432 


0.15000 


0.22500 


7 


16.00000 


21.47563 


0.23559 


0.35339 


8 


18.66666 


23.44431 


0.15046 


0.22568 


9 


21.33333 


25.24850 


0. 15008 


0.22512 


10 


23.99998 


26.98657 


0.15068 


0.22603 


11 


26.66666 


28.49719 


0.15040 


0.22561 


12 


29.33331 


30. 15172 


0. 15009 


0.22514 


13 


31.99998 


31.99998 


0.15012 


0.22517 



E) OBJECTIVE FUNCTION: 



TOTAL STRUCTURE VOLUME: 10.8A1561 



NODE NORMAL STRESS 

1 32852.1 

2 32835.8 

3 32805.9 

A 32730.9 

5 32878.3 

6 32878.2 

7 2103A.9 

8 33175.7 

9 33301.1 

10 33157.7 

11 33220.8 

12 33293.2 

13 33272.5 

F) BOUNDARY CONDITIONS: 

NODE X-DISPL 

1 1 

13 1 

G) SOLUTION vector: 

NODE X-DISPL 

1 O.OOOOOOE+00 

2 0.292592E-01 

3 0.52093AE-01 
A 0.676226E-01 

5 0.806200E-01 

6 0.903615E-01 

7 0.912367E-01 



BENDING STRESS 

0.1 

8291.1 

19268.1 
6A63.6 
7A78.3 
91A3. 1 

31050.3 

16387.2 
2650.6 

3655.2 
15609.0 

18737.4 

0 . 1 



SLOPE 

0 

0 



SLOPE 
-0.927029E-02 
-0.844808E-02 
-0.573258E-02 
-0.436322E-02 
-0.426672E-02 
-0.259323E-02 
0. 155546E-02 



Y-DISPL 

1 

1 



Y-DISPL 

O.OOOOOOE+00 

-0.279060E-01 

-0.511700E-01 

-0.677945E-01 

-0.837384E-01 

-0.969766E-01 

-0.102322E+00 



TOTAL 

32852.2 

41126.9 
52074.0 

39194.5 

40356.5 
42021.4 

52085.2 

49562.9 

35951.7 

36812.9 

48829.8 

52030.6 

33272.6 
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8 

9 

10 

11 

12 

13 



0.815972E-01 

0.675751E-01 

0.529933E-01 

0.398219E-01 

0.227118E-01 

O.OOOOOOE+00 



-0.9^0706E-01 

-0.7971A2E-01 

-0.638010E-01 

-0.A7A255E-01 

-0.26^506E-01 

O.OOOOOOE+00 



0.528782E-02 
0.665123E-02 
0.657931E-02 
0.739218E-02 
0.978718E-02 
0. 111381E-01 
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CASE4 



A 



OPTIMIZATION NO. <4 
OPTIMIZATION SOLUTION 



A) PROBLEM PARAMETERS: 

HORIZ. span: 32.000 

VERT. SPAN : 32.000 

NO OF DESIGN VAR: 2^ 

B) DERIVED CONSTANTS: 

NO OF SYSTEM NODAL POINTS.. 
NO OF DEGREES OF FREEDOM... 
HORIZ. LENGTH PER ELEMENT.. 

C) STRUCTURE LOADING: 

FX 

FY 

FM 

FA 

CONCENTRATED LOAD AT NODE... 



YOUNGS modulus: 30000000.0 
YIELD STRENGTH: 52000.0 
NO OF ELEMENTS: 12 



13 

39 

2.6667 



0.0000 

0.0000 

6000.0000 

0.0000 

7 



DJ ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 



NODE 


X-COORD 


Y-COORD 


BASE 


AREA 


1 


0.00000 


0.00000 


0.15087 


0.22631 


2 


2.66667 


2.69333 


0.15000 


0.22500 


3 


5.33333 


5.38662 


0. 15000 


0.22500 


4 


8.00000 


8.07996 


0. 15000 


0.22500 


5 


10.66666 


10.77329 


0. 15000 


0.22500 


6 


13.33333 


13.46662 


0. 15000 


0.22500 


7 


16.00000 


16. 10661 


0.15377 


0.23066 


8 


18.66666 


18.74661 


0.15000 


0.22500 


9 


21.33333 


21.43994 


0. 15014 


0.22521 


10 


23.99998 


24. 13327 


0. 15005 


0.22507 


11 


26.66666 


26.77327 


0. 15000 


0.22500 


12 


29.33331 


29.41327 


0. 15000 


0.22500 


13 


31.99998 


31.99998 


0. 15000 


0.22500 



E) OBJECTIVE FUNCTION: 



TOTAL STRUCTURE VOLUME: 10.207288 



NODE 


NORMAL STRESS BENDING 


STRESS 


TOTAL 


1 


18.9 




0.0 


18.9 


2 


19.0 


8935.8 


8954.8 


3 


19.0 


17871.6 


17890.6 


4 


19.0 


26806.6 


26825.6 


5 


19.1 


35742. 1 


35761.2 


6 


22.0 


44677.8 


44699.8 


7 


24.3 


52024.6 


52048.9 


8 


22.0 


44303.3 


44325.2 


9 


19.0 


35335.9 


35354.9 


10 


21.8 


26425.4 


26447.2 


11 


24.6 


17591. 1 


17615.7 


12 


27.6 


8749.7 


8777.3 


13 


30.6 




0.1 


30.7 


1 BOUNDARY CONDITIONS: 








NODE 


X-DISPL 


Y-DISPL 


SLOPE 




1 


1 


1 


0 




13 


1 


1 


0 




1 SOLUTION vector: 








NODE 


X-DISPL 


Y-DISPL 


SLOPE 




1 


O.OOOOOOE+00 0, 


OOOOOOE+00 


-0.911602E-02 




2 


0.238804E-01 -0, 


236406E-01 


-0.836558E-02 




3 


0.437104E-01 -0, 


,432711E-01 


-0.610774E-02 




4 


0.554325E-01 -0, 


548738E-01 


-0.234471E-02 




5 


0.549926E-01 -0, 


544349E-01 


0.292347E-02 




6 


0.383368E-01 -0. 


379406E-01 


0.969686E-02 




7 


0.238511E-02 -0. 


162140E-02 


0. 177837E-01 
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FILE: FILE 



CASE4 



A 



8 - 
9 - 
10 - 
11 - 
12 - 
13 



0.336A80E-01 

0.505593E-01 

0.51A326E-01 

0.A05827E-01 

0.21988AE-01 

O.OOOOOOE+00 



0.3A7801E-01 

0.515273E-01 

0.523953E-01 

0.A1AA02E-01 

0.226625E-01 

O.OOOOOOE-i-00 



0.975857E-02 

0.305128E-02 

-0.21508AE-02 

-0.582135E-02 

-0.801783E-02 

-0.87A020E-02 
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file: file 



CASES 



A 



OPTIMIZATION NO. 5 
OPTIMIZATION SOLUTION 



A) PROBLEM parameters: 

HORIZ. SPAN: 32.000 

VERT. SPAN : 32.000 

NO OF DESIGN VAR: 2A 

B) DERIVED constants: 

NO OF SYSTEM NODAL POINTS.. 
NO OF DEGREES OF FREEDOM... 
HORIZ. LENGTH PER ELEMENT.. 

C) STRUCTURE LOADING: 

FX 

FY 

FM 

FA 

CONCENTRATED LOAD AT NODE... 



YOUNGS MODULUS: 30000000.0 
YIELD STRENGTH: 52000.0 
NO OF ELEMENTS: 12 



13 

39 

2.6667 



0.0000 

0.0000 

0.0000 

- 100.0000 

7 



D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 



NODE 


X-COORD 


Y-COORD 


BASE 


AREA 


1 


0.00000 


0.00000 


0.15000 


0.22500 


2 


2.66667 


10.66666 


0. 167A6 


0.25118 


3 


5.33333 


16.A3103 


0.15000 


0.22500 


A 


8.00000 


20.18A51 


0. 15000 


0.22500 


5 


10.66666 


22.93553 


0. 15000 


0.22500 


6 


13.33333 


25. 15350 


0. 15006 


0.22510 


7 


16.00000 


26.92931 


0. 15000 


0.22500 


8 


18.66666 


28.3AA77 


0.15000 


0.22500 


9 


21.33333 


29.A7061 


0. 15000 


0.22500 


10 


23.99998 


30.29001 


0. 15000 


0.22500 


11 


26.66666 


30.90901 


0. 15000 


0.22500 


12 


29.33331 


31.A3A95 


0. 15000 


0.22500 


13 


31.99998 


31.99998 


0. 15000 


0.22500 



E) OBJECTIVE function: 



TOTAL structure VOLUME: 11.32869A 



NODE 


NORMAL STRESS 


BENDING STRESS 


TOTAL 


1 


13798.3 


0.1 


13798.5 


2 


12371.2 


396A1.7 


52012.9 


3 


13857.5 


38158. A 


52015.9 


A 


13897.5 


36095.7 


A9993.2 


5 


13901.5 


38115.7 


52017.2 


6 


13895.2 


38108.5 


52003.6 


7 


13903.2 


3809A.6 


51997.8 


8 


13906.0 


38087.9 


51993.8 


9 


13909. A 


37119. A 


51028.8 


10 


13912.9 


37831.8 


517AA.7 


11 


1392A.A 


35752.2 


A9676.6 


12 


13922.3 


25555.1 


39A77.A 


13 


13910.7 


1.5 


13912.2 



F) BOUNDARY CONDITIONS: 
NODE X-DISPL 


Y-DISPL 


1 


0 


1 


13 


1 


1 



SLOPE 

0 

0 



G) SOLUTION vector: 
NODE X-DISPL 

1 -0.59S168E+00 

2 -0.17A663E+00 

3 -0.196828E-01 
A 0.A62A21E-01 

5 0.752630E-01 

6 0.8A83S8E-01 

7 0.822190E-01 



Y-DISPL 
0.000000E4-00 
-0. 110052E4-00 
-0.184795E+00 
-0.23A2A6E+00 
-0.26A8S0E+00 
-0.278872E+00 
-0.277620E+00 



SLOPE 

-0.A29371E-01 

-0.327189E-01 

-0.217266E-01 

-0.1A1289E-01 

-0.7810A0E-02 

-0.193522E-02 

0.3A902AE-02 
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file: file 



CASES 



A 



8 0.72A252E-01 -0. 26215AE+00 

9 0.587701E-01 -0. 233261E+00 

10 0.AA62A1E-01 -0. 191628E+00 

11 0.30793^E-01 -0. 137662E+00 

12 0.166607E-01 -0. 725301E-01 

13 O.OOOOOOE+OO O.OOOOOOE^fOO 



0.8601A2E-02 

0.13A390E-01 

0.180855E-01 

0.225620E-01 

0.262650E-01 

0.278131E-01 
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file: file 



CASE6 



A 



OPTIMIZATION NO. 6 
OPTIMIZATION SOLUTION 



A) PROBLEM PARAMETERS: 

HORIZ. SPAN; 32.000 
VERT. SPAN : 32.000 
NO OF DESIGN VAR: 2A 



YOUNGS MODULUS: 30000000.0 
YIELD STRENGTH: 52000.0 
NO OF elements: 12 



B) DERIVED constants: 

NO OF SYSTEM NODAL POINTS... 13 

NO OF DEGREES OF FREEDOM.... 39 

HORIZ. LENGTH PER ELEMENT... 2.6667 



C) STRUCTURE LOADING*. 

FX 0. 0000 

FY 0.0000 

FM 0.0000 

FA -100.0000 

CONCENTRATED LOAD AT NODE... 7 



D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 



NODE 


X-COORD 


Y-COORD 


BASE 


AREA 


1 


0.00000 


0.00000 


0. 15000 


0.22500 


2 


2.66667 


6.78236 


0.30329 


0.A5A9A 


3 


5.33333 


lA. 31105 


0.2A58A 


0.36876 


A 


8.00000 


19.67329 


0.15000 


0.22500 


5 


10.66666 


22.9A662 


0. 15000 


0.22500 


6 


13.33333 


2A. 93712 


0. 15000 


0.22500 


7 


16.00000 


26.01016 


0. 15000 


0.22500 


8 


18.66666 


27. 1A098 


0. 15000 


0.22500 


9 


21.33333 


28.00891 


0.15000 


0.22500 


10 


23.99998 


28.76320 


0.15000 


0.22500 


11 


26.66666 


29.500A7 


0.15000 


0.22500 


12 


29.33331 


30.3016A 


0. 1500A 


0. 22506 


13 


31.99998 


31.99998 


0.A3A13 


0.65119 



E) OBJECTIVE FUNCTION: 



TOTAL STRUCTURE VOLUME: 1A.A30920 



NODE 


NORMAL STRESS 


BENDING STRESS 


TOTAL 


1 


1222A.7 


0.1 


1222A.8 


2 


6057. A 


A59A5.8 


52003. 1 


3 


75A7.9 


AAA69.5 


52017. A 


A 


12709.8 


19397.3 


32107. 1 


5 


13071.6 


8753. 1 


2182A.7 


6 


13083. A 


10261.1 


233AA.6 


7 


12983. A 


1171A.5 


2A698.0 


8 


12980.5 


16215.1 


29195.6 


9 


12969.5 


18953.1 


31922.6 


10 


12975.1 


13532.8 


26507.9 


11 


12951.8 


A652.1 


17603.9 


12 


1216A.7 


39871. A 


52036.1 


13 


39A2.6 


A8101.8 


520AA.A 



F) BOUNDARY CONDITIONS: 
NODE X-DISPL 


Y-DISPL 


1 


0 


1 


13 


1 


1 



SLOPE 

0 

1 



G) SOLUTION vector: 
NODE X-DISPL 

1 -0.238121E+00 

2 -0.A19812E-01 

3 0.630188E-01 
A 0.662A82E-01 

5 0.522159E-01 

6 0.AAA5A8E-01 

7 0.A07255E-01 



Y-DISPL 
O.OOOOOOE+00 
-0.792292E-01 
-0.11831AE+00 
-0.122027E+00 
-0. 1129A6E+00 
-0.10A995E+00 
-0.990568E-01 



SLOPE 

-0.323AA2E-01 

-0.223869E-01 

-0.631173E-02 

0.299570E-02 

0.399A36E-02 

0.258829E-02 

0.26811AE-02 
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FILE: FILE 



CASE6 



A 



8 0.355767E-01 -0. 901286E-01 

9 0.296101E-01 -0.757135E-01 

10 0.226288E-01 -0.55A353E-01 

11 0.1A7A3AE-01 -0.31A062E-01 

12 0.702276E-02 -0. 987708E-02 

13 O.OOOOOOE+00 O.OOOOOOE+00 



0.AA789AE-02 

0.667067E-02 

0.867125E-02 

0.921718E-02 

0.6A6195E-02 

O.OOOOOOE+00 
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FILE: FILE 



CASE? 



A 



OPTIMIZATION NO. 7 
OPTIMIZATION SOLUTION 



A) PROBLEM PARAMETERS: 

HORIZ. SPAN: 32.000 
VERT. SPAN : 32.000 
NO OF DESIGN VAR: 2^ 



YOUNGS MODULUS: 30000000.0 
YIELD STRENGTH: 52000.0 
NO OF ELEMENTS: 12 



B) DERIVED CONSTANTS: 

NO OF SYSTEM NODAL POINTS... 13 

NO OF DEGREES OF FREEDOM 39 

HORIZ. LENGTH PER ELEMENT... 2.6667 



C) STRUCTURE LOADING: 

FX 0.0000 

FY 0.0000 

FM 0.0000 

FA “100.0000 

CONCENTRATED LOAD AT NODE... 7 



D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 



NODE 


X-COORD 


Y-COORD 


EASE 


AREA 


1 


0.00000 


0.00000 


1.31120 


3.93359 


2 


2.66667 


2.66377 


1.10192 


3.30575 


3 


5.33333 


5.25856 


0.91330 


2.73989 


4 


8.00000 


7.97715 


0.73940 


2.21821 


5 


10.66666 


11.05349 


0.57427 


1.72280 


6 


13.33333 


14.02065 


0.43146 


1.29439 


7 


16. 00000 


16.77490 


0.31308 


0.93925 


8 


18.66666 


19.50882 


0.30000 


0.90000 


9 


21.33333 


22. 17821 


0.30000 


0.90000 


10 


23.99998 


24.47510 


0.30000 


0.90000 


11 


26.66666 


27.00243 


0.30000 


0.90000 


12 


29.33331 


29.40999 


0.30000 


0.90000 


13 


31.99998 


31.99998 


0.30000 


0.90000 



E) OBJECTIVE function: 



TOTAL STRUCTURE VOLUME: 72.614471 



NODE 


NORMAL STRESS 


1 


0.5 


2 


8.9 


3 


16. 1 


4 


61.1 


5 


123.9 


6 


107.4 


7 


92.3 


8 


79.6 


9 


60.3 


10 


28.9 


11 


9.9 


12 


7.7 


13 


0.0 


F) BOUNDARY CONDITIONS: 


NODE 


X-DISPL Y-: 


1 


1 


13 


0 



BENDING STRESS 


TOTAL 


52070.8 


52071.3 


52066.6 


52075.5 


52056.2 


52072.2 


51982.5 


52043.6 


51879.7 


52003.6 


51886.1 


51993.4 


51919.3 


52011.6 


37074.2 


37153.8 


23343.8 


23404.1 


13394.2 


13423.0 


5932.8 


5942.7 


1530.9 


1538.6 


5.7 


5.7 



•L SLOPE 

1 1 

0 0 



G) SOLUTION VECTOR: 
NODE X-DISPL 

1 O.OOOOOOE+00 

2 0.597678E-02 

3 0.230543E-01 

4 0.528078E-01 

5 0.100533E+00 

6 0.160387E+00 

7 0.228394E+00 



Y-DISPL 

O.OOOOOOE+00 

-0.598316E-02 

-0.235304E-01 

-0.527179E-01 

-0.941091E-01 

-0.147924E+00 

-0.213783E+00 



SLOPE 

O.OOOOOOE+00 

-0.436127E-02 

-0.866592E-02 

-0.130684E-01 

-0.177673E-01 

-0.223668E-01 

“0.267883E-01 
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FILE: FILE 



CASE? 



A 



8 

9 

10 

11 

12 

13 



0.307128E+00 
0.39238AE+00 
0. A702A1E+00 
0.558680E4-00 
0.6AA186E+00 
0.736590E+00 



-0,290595E+00 

-0.375776E+00 

-0.A66157E+00 

-0.559A72E+00 

-0.65A177E+00 

-0.7A9316E+00 



-0.305781E-01 

-0.331111E-01 

-0.3«A77E-01 

-0.353367E-01 

-0.3563^9E-01 

-0.356983E-01 
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FILE: FILE 



CASES 



A 



OPTIMIZATION NO. 8 
OPTIMIZATION SOLUTION 



A) PROBLEM PARAMETERS: 

HORIZ. SPAN: 32.000 

VERT. SPAN : 32.000 

NO OF DESIGN VAR: 2A 

B) DERIVED CONSTANTS; 

NO OF SYSTEM NODAL POINTS.. 
NO OF DEGREES OF FREEDOM... 
HORIZ. LENGTH PER ELEMENT.. 

C) STRUCTURE LOADING: 

FX 

FY 

FM 

FA 

CONCENTRATED LOAD AT NODE... 



YOUNGS MODULUS: 30000000.0 
YIELD strength: 52000.0 
NO OF ELEMENTS: 12 



13 

39 

2.6667 



0.0000 

0.0000 

0.0000 

- 100.0000 

7 



D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 



NODE 


X-COORD 


Y-COORD 


EASE 


AREA 


1 


0.00000 


0.00000 


0.A029A 


0.60AA0 


2 


2.66667 


3.27370 


0.15000 


0.22500 


3 


5.33333 


8.18358 


0.15000 


0.22500 


A 


8.00000 


12.A9695 


0.15000 


0.22500 


5 


10.66666 


15.68683 


0.15000 


0.22500 


6 


13.33333 


18.A2828 


0.15000 


0.22500 


7 


16.00000 


20.901A0 


0. 15000 


0.22500 


8 


18.66666 


23.20117 


0.15000 


0.22500 


9 


21.33333 


25.28195 


0.15000 


0.22500 


10 


23.99998 


27.28166 


0.15000 


0.22500 


11 


26.66666 


29.15797 


0.15000 


0.22500 


12 


29.33331 


30.5A808 


0.15000 


0.22500 


13 


31.99998 


31.99998 


0.21518 


0.32277 



E) OBJECTIVE FUNCTION: 



TOTAL STRUCTURE VOLUME: 11.308629 



NODE NORMAL STRESS 

1 9208.7 

2 2508^.6 

3 25^»A3. 1 

A 25A60.7 

5 25A36.8 

6 25386.7 

7 25363. A 

8 25365.1 

9 25379.9 

10 25A1A.8 

11 2551A.8 

12 25579.6 

13 17825.3 

F) BOUNDARY CONDITIONS: 

NODE X-DISPL 

1 1 

13 1 

G) SOLUTION vector: 

NODE X-DISPL 

1 O.OOOOOOE+00 

2 0.113937E-01 

3 0.3A2988E-01 
A 0.5051A1E-01 

5 0.6A8279E-01 

6 0.7A9570E-01 

7 0.768717E-01 



BENDING STRESS 
A2801.0 
3375. A 
6676.3 
9375.9 
3831.6 
181A7.5 
266A9.6 
25377.8 

20239.0 
1767.0 

25373. A 

26397.0 
3AA01.1 



SLOPE 

1 

1 



SLOPE 
O.OOOOOOE+00 
-0.568126E-02 
-0.AA3323E-02 
-0.A737AAE-02 
S2A971E-02 
-0.338176E-02 
0.2388A5E-03 



Y-DISPL 

1 

1 



Y-DISPL 
0.000000E-l>00 
-0.117173E-01 
-0.295A78E-01 
-0.AA6311E-01 
-0.611976E-01 -0. 
-0.755686E-01 
-0.821558E-01 



TOTAL 

52009.6 
28A60.0 

32119.3 
3A836.6 

29268.3 
A353A.2 
52012.9 
507A2.9 
A5618.8 
27181.8 
50888.2 

51976.7 
52226. A 



1 01 



file: file 



CASES 



A 



8 0.693A85E-01 -0. 779902E-01 

9 0.5AA237E-01 -0. 635129E-01 

10 0.3A6055E-01 -0. A17868E-01 

11 0.1575A9E-01 -0. 198010E-01 

12 0.522320E-02 -0. 51A698E-02 

13 O.OOOOOOE+00 O.OOOOOOE+00 



0.A31015E-02 

0.773890E-02 

0.936886E-02 

0.765835E-02 

0.A19866E-02 

O.OOOOOOE+00 



1 02 



file: file 



CASE9 



A 



OPTIMIZATION NO. 9 
OPTIMIZATION SOLUTION 



A) PROBLEM parameters: 

HORIZ. SPAN: 32,000 

VERT. SPAN ; 18.A75 

NO OF DESIGN VAR: 2A 

B) DERIVED CONSTANTS: 

NO OF SYSTEM NODAL POINTS.. 
NO OF DEGREES OF FREEDOM... 
HORIZ. LENGTH PER ELEMENT.. 

C) STRUCTURE LOADING: 

FX 

FY 

FM 

FA 

CONCENTRATED LOAD AT NODE... 



YOUNGS MODULUS: 30000000.0 
YIELD STRENGTH: 52000.0 
NO OF elements: 12 



13 

39 

2.6667 



0.0000 

0.0000 

0.0000 

- 100.0000 

7 



D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 



NODE 


X-COORD 


Y-COORD 


BASE 


AREA 


1 


0.00000 


0.00000 


0. 15003 


0.2250A 


2 


2.66667 


2.28090 


0.15000 


0.22500 


3 


5.33333 


A. 52538 


0.15025 


0.22537 


A 


8.00000 


6.5736A 


0.15012 


0.22518 


5 


10.66666 


8.AA217 


0.15017 


0.22526 


6 


13.33333 


10.1A7A6 


0. 150 lA 


0.22521 


7 


16.00000 


11.701A6 


0.15023 


0.22535 


8 


18.66666 


13.11621 


0.15000 


0.22500 


9 


21.33333 


1A.397AA 


0.15005 


0.22507 


10 


23.99998 


15.55AA5 


0.1502A 


0.22535 


11 


26.66666 


16.59A73 


0.15000 


0.22500 


12 


29.33331 


17.5263A 


0.15000 


0.22500 


13 


31.99998 


18.A7A99 


0. 15000 


0.22500 



E) OBJECTIVE FUNCTION: 



TOTAL STRUCTURE VOLUME: 8.390676 



NODE 


NORMAL STRESS 


BENDING STRESS 


TOTAL 


1 


33701.0 


0.1 


33701.1 


2 


3370A.7 


182A2.0 


519A6.6 


3 


336A6.0 


18393.0 


52039.0 


A 


33675.5 


18362.1 


52037.6 


5 


33663.0 


18380.3 


520A3.3 


6 


33671.2 


18369.3 


520A0.5 


7 


336A9.9 


18380.3 


52030.2 


8 


33702.3 


18278.8 


51981.2 


9 


33692.0 


18333.3 


52025.3 


10 


336A9.5 


18317.5 


51967.0 


11 


33703. A 


18112.9 


51816.3 


12 


33701.9 


17257.1 


50959.1 


13 


33699.5 


0.1 


33699.6 



F) BOUNDARY CONDITIONS: 
NODE X-DISPL 


Y-DISPL 


1 


1 


1 


13 


1 


1 



SLOPE 

0 

0 



G) SOLUTION VECTOR: 
NODE X-DISPL 

1 O.OOOOOOE+00 

2 0.273039E-01 

3 0.A88191E-01 
A 0.62A682E-01 

5 0.696080E-01 

6 0.713910E-01 

7 0.6878A6E-01 



Y-DISPL 

O.OOOOOOE+OO 

-0.379869E-01 

-0.6962A8E-01 

-0.93588AE-01 

-0.1101A6E+00 

-0.119528E+00 

-0.121933E+00 



SLOPE 

-0.137583E-01 

-0.123359E-01 

-0.9A983AE-02 

-0.67518AE-02 

-0.A0931AE-02 

-0.150820E-02 

0.1012A1E-02 
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FILE: FILE 



CASE9 



A 



8 0.626172E-01 -0. 1 17539E+00 

9 0.536328E-01 -0. 106513E+00 

10 0.A2A687E-01 -0. 889785E-01 

11 0.296891E-01 -0. 650603E-01 

12 0.15809AE-01 -0. 3A9536E-01 

13 O.OOOOOOE+00 O.OOOOOOE't'OO 



0.3A7161E-02 
0.587863E-02 
0.82^i613E-02 
0. 10563AE-01 
0. 127837E-01 
0. 138690E-01 



1 04 



FILE: FILE 



CASEIO A 



OPTIMIZATION NO. 10 
OPTIMIZATION SOLUTION 



A) PROBLEM PARAMETERS: 

HORIZ. SPAN: 32.000 
VERT. SPAN : 55.^26 
NO OF DESIGN VAR: 2A 



YOUNGS modulus: 30000000.0 
YIELD STRENGTH: 52000.0 
NO OF elements: 12 



B) DERIVED constants: 

NO OF SYSTEM NODAL POINTS. . . 13 

NO OF DEGREES OF FREEDOM 39 

HORIZ. LENGTH PER ELEMENT... 2.6667 



C) STRUCTURE LOADING: 

FX 0.0000 

FY 0.0000 

FM 0.0000 

FA -100.0000 

CONCENTRATED LOAD AT NODE... 7 



D) ELEMENTAL DIMENSIONS AND STRESS DISTRIBUTION: 



NODE 


X-COORD 


Y-COORD 


BASE 


AREA 


1 


0.00000 


0.00000 


0.15000 


0.22500 


2 


2.66667 


11.92522 


0.15000 


0.22500 


3 


5.33333 


19.67612 


0. 15000 


0.22500 


A 


8.00000 


26.02310 


0. 15000 


0.22500 


5 


10.66666 


30.8803A 


0.15000 


0.22500 


6 


13.33333 


3A. 78506 


0. 20A05 


0.30607 


7 


16. 00000 


39. 13298 


0.15000 


0.22500 


8 


18.66666 


A2.259A3 


0.15000 


0.22500 


9 


21.33333 


A5. 35509 


0. 15000 


0.22500 


10 


23.99998 


A8. 1760A 


0.15000 


0.22500 


11 


26.66666 


50.91873 


0.15000 


0.22500 


12 


29.33331 


53.12798 


0.15000 


0.22500 


13 


31.99998 


55.A2598 


0. 15000 


0.22500 



E) OBJECTIVE FUNCTION: 



NODE 


TOTAL STRUCTURE VOLUME: 15. 1A862A 

NORMAL STRESS BENDING STRESS 


TOTAL 


1 


37992.2 


0. 1 


37992.3 


2 


38009.5 


13993.7 


52003.2 


3 


380A5.6 


13552.0 


51597.5 


A 


3808A.5 


9528.0 


A7612.5 


5 


38055.6 


1160.7 


39216.3 


6 


27936.6 


2A086.0 


52022.6 


7 


38030.6 


1322A.5 


51255.2 


8 


38061. A 


13822.2 


51883.6 


9 


38059.8 


13996.9 


52056.7 


10 


3806A.1 


1397A. 1 


52038. 1 


11 


38067. A 


A071.1 


A2138.A 


12 


38068.1 


139A6.2 


5201A.3 


13 


38069.7 


0. 1 


38069.9 


F) BOUNDARY CONDITIONS: 

NODE X-DISPL Y-DISPL 


SLOPE 




1 


1 1 


0 




13 


1 1 


0 





G) SOLUTION VECTOR; 
NODE X-DISPL 

1 O.OOOOOOE+00 

2 0.103071E+00 

3 0.129693E+00 
A 0.129351E+00 

5 0.133686E+00 

6 0.135^01E+00 

7 0.122050E+00 



Y-DISPL 

O.OOOOOOE+00 

-0.389057E-01 

-0.590526E-01 

-0.683836E-01 

-0.787921E-01 

-0.861105E-01 

-0.8A3AA0E-01 



SLOPE 

-0.101930E-01 

-0.639297E-02 

-0.1375A8E-02 

-0.759869E-03 

-0.207601E-02 

0.737877E-03 

0.261^56E-02 
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8 0.11169AE+00 -0.823635E-01 

9 0.961<»71E-01 -0.758125E-01 

10 0.7A70A9E-01 -0« 623201E-01 

n 0.A85993E-01 -0. A37086E-01 

12 0.262193E-01 -0. 23582AE-01 

13 O.OOOOOOE+00 O.OOOOOOE+00 



0.266911E-02 
0.519A96E-02 
0.760786E-02 
0.8AA971E-02 
0.920963E-02 
0. 103006E-01 
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